Harmonic  Approaches  to  Non-Intrusive  Load  Diagnostics 

by 

Ashley  E.  Fuller 

B.S.,  Electrical  Engineering  &  Technology  (2003),  Old  Dominion  University 

Submitted  to  the  Department  of  Mechanical  Engineering  in  Partial  Fulfillment  of  the 

Requirements  for  the  Degrees  of 

Naval  Engineer 
and 

Master  of  Science„imSystejn_Design  and  Management 

at  the 

Massachusetts  Institute  of  Technology,  June  2008 


©  2008  Ashley  E.  Fuller.  All  rights  reserved. 


The  author  hereby  grants  to  MIT  permission  to  reproduce  and  to  distribute  publicly  paper 
and  electronic  copies  of  this  thesis  document  in  whole  or  in  part  in  any  medium  now 

knovm  or  hereafter  created. 


Signature  of  Author_ 


Certified  by 


(MyJ-  Qm 


Department  of  Mechanical  Engineering 

May  9,  2008 


Robert  W.  Cox,  Assistant  Professor  University  of  North  Carolina  Charlotte 

Department  of  Electrical  and  Computer  Engineering 
^0-  Thesis  Advisor 

Certified  by  ^ _ _ _ 

Steven  B.  Leeb 

Professor  of  Electrical  Engineering  and  Computer  Science  &  Mechanical  Engineering 
Departments  of  Elps^rical  Engineering  ^d^omputer  Science  &  Mechanical  Engineering 

Thesis  Advisor 

4  1/  V  \  /  \/ 

Accepted  by 


Pat  Hale,  Dh 


Accepted  by 


ystem  Design  and  Management  Fellows  Program 
Engineering  Systems  Division 
Thesis  Advisor 


Lallit  AnanoTPf^^ssor  of  Mechanical  Engineering 
Chairman,  Department  Committee  on  Graduate  Students 
Department  of  Mechanical  Engineering 


1 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports 
(0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be 
subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 


1.  REPOKT  DAli.  (DD-MM-YYYY) 

XX-06-2008 


4.  TITLE  AND  SUBTITLE 

Harmonic  Approaches  to  Non-Intrusive  Load  Diagnostics 


Master's  Thesis 


3.  DATES  COVERED  (From  -  To) 

JAN-JUN  2008 


5a.  CONTRACT  NUMBER 

N62271-97-G-0026 

5b.  GRANT  NUMBER 


PROGRAM  ELEMENT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Massachusetts  Institute  of  Technology 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 
Monterey,  CA  93943 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 


1.  DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


14.  ABSTRACT 

The  Non-Intrusive  Load  Monitor  (NILM)  is  a  system  that  monitors,  records  and  processes 
voltage  and  current  measurements  to  establish  the  operating  characteristics  of  individual  loads  on 
a  load  center  from  a  single  aggregate  measurement.  The  NILM  can  also  be  used  to  actively 
monitor  degradation  or  diagnose  specific  system  failures.  Current  NILM  research  conducted  at 
the  Massachusetts  Institute  of  Technology’s  Laboratory  for  Electromagnetic  and  Electronic 
Systems  (LEES)  is  exploring  the  application  and  expansion  of  NILM  technology  for  the  use  of 
monitoring  a  myriad  of  electromechanical  loads.  This  thesis  presents  a  fundamental  guide  to 
understanding  NILM  operation  using  laboratory  bench  testing  and  demonstrates  its  potential  to 
detect  an  array  of  electric  machine  failures  before  they  become  catastrophic. 


16.  SECURITY  CLASSIFICATION  OF: 
a.  REPORT  I  b.  ABSTRACT  I  c.  THIS  PAGE 


17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 
PAGES  Tibbitts,  Educational  Technician 

19b.  TELEPHONE  NUMBER  (Include  area  code) 
(831)  656-2319  civins@nps.edu 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


INSTRUCTIONS  FOR  COMPLETING  SF  298 


(1  ■  REPORT  DATEj  Full  publication  date,  including 
day,  month,  if  available.  Must  cite  at  least  the  year 
and  be  Year  2000  compliant,  e.g.  30-06-1998; 
xx-06-1998;  xx-xx-1998. 

(2.  REPORT  TYP^  State  the  type  of  report,  such  as 
final,  technical,  interim,  memorandum,  master's 
thesis,  progress,  quarterly,  research,  special,  group 
study,  etc. 

(3.  DATES  COVERED)  Indicate  the  time  during 
which  the  work  was  performed  and  the  report  was 
written,  e.g.,  Jun  1997  -  Jun  1998;  1-10  Jun  1996; 
May  -  Nov  1998;  Nov  1998. 

(4.  TITLE}  Enter  title  and  subtitle  with  volume 
number  and  part  number,  if  applicable.  On  classified 
documents,  enter  the  title  classification  in 
parentheses. 

(5a.  CONTRACT  NUMBER  Enter  all  contract 
numbers  as  they  appear  in  the  report,  e.g. 
F33615-86-C-5169. 

5b.  GRANT  NUMBER.  Enter  all  grant  numbers  as 
they  appear  in  the  report,  e.g.  AFOSR-82-1 234. 

5c.  PROGRAM  ELEMENT  NUMBER.  Enter  all 
program  element  numbers  as  they  appear  in  the 
report,  e.g.  61 101  A. 

5d.  PROJECT  NUMBER.  Enter  all  project  numbers 
as  they  appear  in  the  report,  e.g.  1 F665702D1 257; 
ILIR. 

5e.  TASK  NUMBER.  Enter  all  task  numbers  as  they 
appear  in  the  report,  e.g.  05;  RF0330201;  T4112. 

5f.  WORK  UNIT  NUMBER.  Enter  all  work  unit 
numbers  as  they  appear  in  the  report,  e.g.  001; 
AFAPL30480105. 


(6.  AUTHOR(S))  Enter  name(s)  of  person(s) 
responsible  for  writing  the  report,  performing  the 
research,  or  credited  with  the  content  of  the  report. 
The  form  of  entry  is  the  last  name,  first  name,  middle 
initial,  and  additional  qualifiers  separated  by  commas, 
e.g.  Smith,  Richard,  J,  Jr. 

(7.  PERFORMING  ORGANIZATION  NAME(S)  AND) 

(ADDRESS(ES).)  Self-explanatory. 


8.  PERFORMING  ORGANIZATION  REPORT  NUMBER. 

Enter  all  unique  alphanumeric  report  numbers  assigned 
by  the  performing  organization,  e.g.  BRL-1234; 
AFWL-TR-85-4017-Vol-21-PT-2. 

(9.  SPONSORING/MONITORING  AGENCY  NAME(^ 

(AND  ADDRESS(ES).  )Enter  the  name  and  address  of  the 
organization(s)  financially  responsible  for  and  monitoring 
the  work. 

10.  SPONSOR/MONITOR’S  ACRONYM(S).  Enter,  if 
available,  e.g.  BRL,  ARDEC,  NADC. 

1 1 .  SPONSOR/MONITOR'S  REPORT  NUMBER(S). 

Enter  report  number  as  assigned  by  the  sponsoring/ 
monitoring  agency,  if  available,  e.g.  BRL-TR-829;  -215. 

(12.  DISTRIBUTION/AVAILABILITY  STATEMENT^  Use 

agency-mandated  availability  statements  to  indicate  the 
public  availability  or  distribution  limitations  of  the 
report.  If  additional  limitations/  restrictions  or  special 
markings  are  indicated,  follow  agency  authorization 
procedures,  e.g.  RD/FRD,  PROPIN,  ITAR,  etc.  Include 
copyright  information. 

13.  SUPPLEMENTARY  NOTES.  Enter  information  not 
included  elsewhere  such  as:  prepared  in  cooperation 
with;  translation  of;  report  supersedes;  old  edition 
number,  etc. 

(14.  ABSTRACT)  A  brief  (approximately  200  words) 
factual  summary  of  the  most  significant  information. 

15.  SUBJECT  TERMS.  Key  words  or  phrases 
identifying  major  concepts  in  the  report. 

16.  SECURITY  CLASSIFICATION.  Enter  security 
classification  in  accordance  with  security  classification 
regulations,  e.g.  U,  C,  S,  etc.  If  this  form  contains 
classified  information,  stamp  classification  level  on  the 
top  and  bottom  of  this  page. 

(17.  LIMITATION  OF  ABSTRACT.)  This  block  must  be 
completed  to  assign  a  distribution  limitation  to  the 
abstract.  Enter  UU  (Unclassified  Unlimited)  or  SAR 
(Same  as  Report).  An  entry  in  this  block  is  necessary  if 
the  abstract  is  to  be  limited. 


Standard  Form  298  Back  (Rev.  8/98) 


Page  Intentionally  Left  Blank 


2 


Harmonic  Approaches  to  Non-Intrusive  Load  Diagnostics 

by 

Ashley  E.  Fuller 


Submitted  to  the  Department  of  Mechanical  Engineering  on  May  9,  2008  in  Partial  Fulfillment 

of  the  Requirements  for  the  Degrees  of 

Naval  Engineer 
and 

Master  of  Science  in  System  Design  and  Management 


Abstract 

The  Non-Intrusive  Load  Monitor  (NILM)  is  a  system  that  monitors,  records  and  processes 
voltage  and  current  measurements  to  establish  the  operating  characteristics  of  individual  loads  on 
a  load  center  from  a  single  aggregate  measurement.  The  NILM  can  also  be  used  to  actively 
monitor  degradation  or  diagnose  specific  system  failures.  Current  NILM  research  conducted  at 
the  Massachusetts  Institute  of  Technology’s  Laboratory  for  Electromagnetic  and  Electronic 
Systems  (LEES)  is  exploring  the  application  and  expansion  of  NILM  technology  for  the  use  of 
monitoring  a  myriad  of  electromechanical  loads.  This  thesis  presents  a  fundamental  guide  to 
understanding  NILM  operation  using  laboratory  bench  testing  and  demonstrates  its  potential  to 
detect  an  array  of  electric  machine  failures  before  they  become  catastrophic.  The  NILM’s  ability 
to  the  monitor  the  current  spectrum  of  electric  machines  can  be  used  to  immediately  diagnose 
multiple  common  system  casualties  and  detect  unusual  system  operation.  Clean  current 
spectrum  regions  can  be  exploited  by  selecting  induction  machine  design  characteristics  that 
result  in  eccentric  modulation  frequencies  occurring  in  areas  free  of  supply  frequency  harmonics. 

Current  spectrum  analysis  was  used  to  demonstrate  the  NILM’s  potential  to  monitor  multiple 
machines  from  an  aggregate  source  and  discuss  intersystem  impedances.  It  can  be  shown  that 
multiple  machines  with  slightly  varied  physical  characteristics,  such  as  induction  motor  rotor 
slots,  coupled  with  using  clean  current  spectral  regions  support  automated  diagnostic  system 
development.  Measurements  and  experimentation  were  conducted  in  the  LEES  laboratory  and 
the  Industrial  Support  Center  electric  shop,  Boston. 

Thesis  Advisor:  Steven  B.  Leeb 
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Thesis  Advisor:  Robert  W.  Cox 
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Thesis  Advisor:  Pat  Hale 
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1  Introduction 


1.1  Background 

In  theory,  the  Non-Intrusive  Load  Monitor  or  NILM  can  determine  the  operating 
schedule  of  electrical  loads  in  a  target  system  from  measurements  made  at  a  centralized  location, 
such  as  the  electric  utility  service  entry  [1].  It  has  shown  the  ability  to  identify  the  operation  of 
specific  electromechanical  devices  in  an  expanding  collection  of  distribution  networks.  This 
thesis  establishes  a  bench  test  to  observe  NILM  operation  using  standard  test  equipment. 
Voltage  and  Current  are  simulated  in  order  to  provide  a  known  NILM  input  and  to  confirm  the 
output  at  select  stages  of  NILM  data  processing. 

An  induction  machine  is  used  to  test  the  NILM’s  ability  to  observe  common  induction 
machine  problems.  A  framework  similar  to  that  used  in  vibration  analysis  is  employed  to  locate 
and  explain  physical  eccentricities  associated  with  common  induction  machine  faults.  Initially, 
both  the  vibration  spectrum  and  current  spectrum  are  presented  in  parallel  to  validate  results. 
Realistic  faults  are  inserted  into  a  test  motor,  then  current  and  vibration  harmonics  associated 
with  those  faults  are  confirmed.  This  thesis  identifies  specific  current  spectrum  harmonics  and 
uses  them  to  present  a  potential  diagnostic  framework. 

Current  harmonics  associated  with  the  physical  design  of  an  induction  machine  present 
the  potential  to  differentiate  multiple  motors  from  an  aggregate  source.  Multiple  motors  were 
monitored  using  a  single  sensor.  A  fault  was  inserted  into  one  machine  with  the  goal  of  defining 
a  diagnostic  scheme  that  could  determine  the  faulty  machine  from  an  aggregate  current 
measurement.  The  single  machine  diagnostic  framework  is  then  expanded  to  include  a  multiple 
motor  framework. 

Lastly,  the  NILM  architecture  is  analyzed  using  techniques  introduced  in  System  Design 
and  Management  courses.  A  new  architecture  is  presented  using  fewer  sub-assemblies  and 
increased  capability.  The  new  NILM  architecture  is  constructed,  installed  aboard  a  United  States 
Coast  Guard  Cutter  and  tested  with  an  improved  software  package. 


1.2  NILM  Overview 

Figure  1-1  shows  the  block  diagram  of  a  standard  NILM.  Note  that  the  NILM  measures 
the  aggregate  current  flowing  to  a  bank  of  electrical  loads.  It  then  disaggregates  the  operating 
schedule  of  individual  loads  using  signal-processing  techniques  [2].  In  an  engineering  plant,  the 
candidate  installation  locations  include  generator  output  busses  and  distribution  panels. 
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Figure  1-1:  Diagram  showing  the  fundamental  signal  flow  path  in  a  NILM 


Using  measurements  of  the  line  voltage  and  aggregate  current,  a  software -based 
preprocessor  onboard  the  NILM  computes  time- varying  estimates  of  the  frequency  content  of  the 
measured  line  current  [3].  Formally,  these  time-varying  estimates,  or  spectral  envelopes,  are 
defined  as  the  quantities  [4]; 


m  T  ^ ' 


i(r)sin(mtyr)b?i 


b  u)=  —  ff  ^nr)cos(mtyrVr 

^  V  /  rj.  if- 1  \  /  \  r 


(1-1) 


(1-2) 


These  equations  are  essentially  time  varying  Fourier-series  analysis  equations  evaluated 
over  a  moving  window  of  length  T  [5].  The  coefficients  am{t)  and  bm{t)  contain  time-local 
information  about  the  frequency  content  of  i(t).  Provided  that  the  basis  terms  sin(mty  t)  and 
cos{mcot)  are  synchronized  to  the  line  voltage,  the  spectral  envelope  coefficients  have  a  useful 
physical  interpretation  as  real,  reactive,  and  harmonic  content  [6]. 

The  spectral  envelopes  computed  by  the  preprocessor  are  passed  to  an  event  detector  that 
identifies  the  operation  of  each  of  the  major  loads  on  the  monitored  electrical  service.  In  a 
modern  NILM,  identification  is  performed  using  both  transient  and  steady-state  information  [7], 
Field  studies  have  demonstrated  that  transient  details  are  particularly  powerful  because  the 
transient  electrical  behavior  of  a  particular  load  is  strongly  influenced  by  the  physical  task  that  is 
performed  [8].  As  shown  in  Fig.  1-2,  the  physical  differences  between  an  incandescent  lamp  and 
an  induction  machine  result  in  different  transient  patterns.  Figure  1-3  demonstrates  the  positive 
identification  of  an  induction  motor  driving  a  small  vacuum  pump.  Further  details  of  the 
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detection  and  identification  process  can  be  found  in  Lee  [7]  and  Leeb  [6]. 


I 

< 


0  0.05  0.1  0.15  0.2  0.25 


Time;  Sec 


Figure  1-2:  Top  trace:  Current  drawn  dnring  the  start  of  an  incandescent  lamp.  Bottom  trace:  Stator 
enrrent  drawn  during  the  start  of  an  nnloaded,  fractional  horsepower  induction  machine  [2]. 
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Fignre  1-3.  Measnred  enrrent  and  compnted  power  dnring  the  start  of  1.7hp  vacuum  pump  motor.  Also 
shown  in  the  power  plot  is  a  section  of  the  template  that  has  been  successfully  matched  to  the  observed 
transient  behavior  by  the  NILM’s  event  detector  [2]. 


The  final  block  in  Fig.  1-1  is  the  NILM’s  diagnostics  and  systems  management  module. 
This  software  unit  assesses  load  status  using  any  required  combination  of  current  data,  voltage 
data,  spectral  envelopes,  and  load  operating  schedules  [9].  The  successful  application  of  this 
module  has  been  demonstrated  in  numerous  publications  eg.,  [3],  [7],  [8],  [10].  Shipboard 
applications  are  highlighted  in  [8],  [11],  [12],  and  [13]. 
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As  shown  in  Fig.  1-1,  the  NILM  is  designed  to  interact  with  human  or  automated  devices 
in  a  number  of  different  ways.  For  instance,  the  NILM  can  use  its  diagnostic  information  to 
command  certain  loads  to  either  commence  or  cease  operations.  Additionally,  the  NILM  can 
provide  regular  status  reports  to  human  operators.  To  assist  in  future  maintenance  operations,  the 
NILM  stores  all  of  the  relevant  data  streams  (i.e.  currents,  voltages,  operating  schedules,  etc.)  in 
either  a  local  or  remote  database  [14].  The  NILM’s  storage  capabilities  make  it  possible  for  the 
operator  to  perform  historical  data  trending.  Note  that  this  off-line  analysis  can  be  conducted  on 
a  remote  PC  using  convenient  software  packages  such  as  Microsoft  Excel  [14].  The  following 
sections  describe  how  these  capabilities  can  be  used  to  prevent  or  detect  certain  critical  shipboard 
faults. 


1.3  Objectives  and  Motivation 

The  NILM  is  progressing  toward  production  and  commercialization  to  be  used  in  both 
shipboard  and  commercial  applications.  However,  there  are  several  important  system 
confirmation  tests  that  should  be  performed.  The  NILM  must  be  field  tested  in  systems  where 
diagnostics  have  already  been  developed  and  confirmed.  One  example  of  an  existing  system  test 
that  can  be  used  to  confirm  NILM  employment  is  vibration  analysis,  which  is  a  part  of  the  Coast 
Guard’s  Condition  Based  Maintenance  program.  The  United  States  Coast  Guard  is  currently 
using  vibration  analysis  to  monitor  the  health  of  shipboard  systems  to  determine  system 
maintenance  schedules.  There  are  a  number  of  shortfalls  to  this  system.  First,  it  requires 
someone  familiar  with  the  vibration  monitoring  process  and  appropriate  training  to  travel  to  the 
ship,  configure  the  target  system  for  analysis,  run  the  test  (usually  more  than  once),  record  the 
data,  and  return  to  enter  the  data  in  the  appropriate  data  base,  which  at  the  present  moment  is  an 
Excel  spreadsheet. 

The  vision  of  the  NILM  is  that  a  minimal  number  of  sensors  will  eventually  be  used  to 
monitor  multiple  loads.  The  NILM’s  ability  to  monitor  multiple  loads  using  an  aggregate  current 
measurement  must  be  tested.  Current  versions  of  the  NILM  identify  loads  simply  by  looking  at 
transient  patterns  and  steady-state  power  changes.  This  may  not  be  enough  when  looking  at 
multiple  systems.  One  option  that  will  be  investigated  is  the  use  of  system  spectral  information. 
Analogous  to  the  mechanical  vibration  monitoring,  the  electrical  spectral  information  presents  a 
greater  potential  for  multi-system  remote  monitoring  and  overall  cost  savings. 

This  thesis  will  explain  the  basic  operation  of  each  NILM  system  through  a  discussion  of 
the  internal  signal  processing  and  using  laboratory  bench  test  examples.  A  multi-channel 
waveform  generator  and  signal  generator  are  used  to  modulate  a  signal  representing  current. 
That  current  signal  is  fed  into  the  NILM  and  each  step  of  signal  processing  is  illustrated.  It  will 
assist  in  the  field-validation  study  by  comparing  NILM  with  the  current  accepted  practice  of 
vibration  analysis.  Common  areas  of  vibration  analysis  success,  such  as  induction  machine 
imbalance  and  bearing  failure,  will  be  used  to  confirm  NILM  results.  Induction  machine  spectral 
effects  will  be  used  to  take  a  detailed  look  at  machine  health.  This  thesis  will  catalogue  and 
explain  common  expected  spectral  effects  and  use  them  to  show  the  NILM’s  superior  potential. 
A  procedure  for  monitoring  multiple  induction  machine  loads  will  be  explained  and  feasibility 
investigated  using  the  methods  established  in  the  single  machine  diagnostic  section.  The  second 
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half  of  this  thesis  will  analyze  the  current  NILM  architecture,  define  the  next  generation  NILM 
architecture  and  present  this  new  architecture.  The  NILM  system  architecture  will  be  analyzed 
using  0PM  (Object-Process  Methodology),  system  decomposition,  and  subsystem 
interdependency  analysis  tools  like  DSM  (Design  Structure  Matrices)  presented  in  System 
Design  and  Management  curriculum.  The  next  generation  NILM  will  be  constructed  according 
to  the  output  of  previously  mentioned  tools  and  explained.  Lastly,  a  system  (NILM)  to  system 
(host  or  target)  interface  will  be  examined  to  determine  a  better  installation  methodology. 
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2  Using  the  NILM  to  Monitor  Spectral  Quantities 

2. 1  Introduction 


The  NILM  is  a  monitoring  system  for  extracting  useful  information  about  any  target  or 
host  system  that  uses  electromechanical  devices.  It  has  a  low  installation  cost  and  high  reliability 
because  it  uses  a  bare  minimum  of  sensors.  One  important  feature  of  the  NILM  is  its  ability  to 
detect  the  influence  of  mechanical  effects  on  electromagnetic  actuators  like  motors.  Induction- 
motor  stator  currents,  for  instance,  contain  a  wealth  of  information  about  the  motor’s  mechanical 
load  and  its  operational  status  [15],  [16].  Mechanical  non-idealities  such  as  load  imbalances  and 
rotor  eccentricity  cause  periodic  variations  in  the  amplitude  of  the  current.  Such  variations  can  be 
analyzed  either  by  directly  evaluating  the  frequency  spectrum  of  the  current  or  by  examining  the 
spectrum  of  a  demodulated  signal.  Both  analytical  approaches  are  valuable  because  many 
important  failures  can  be  detected  by  analyzing  current  waveforms. 


2.1.1  Current  Spectrum  Content  in  Induction  Machines 

The  dynamics  in  electromechanical  systems  can  often  be  seen  in  the  frequency  spectrum 
of  electrical  current.  For  example,  mechanical  vibrations  in  induction  machines  can  be 
associated  with  variations  in  the  physical  air  gap  of  the  motor  [18].  These  effects  cause  a  slow 
variation  in  the  torque  produced  by  the  machine,  and  thus  induce  a  corresponding  variation  in  the 
amplitude  of  its  current.  In  effect,  the  current  is  said  to  be  amplitude-modulated.  The  detection 
of  such  modulation  can  be  extremely  useful  from  a  diagnostic  perspective. 

Raw  current  data  obtained  by  the  NILM  can  be  analyzed  to  search  for  variations  such  as 
the  amplitude  modulation  described  above.  To  understand  the  process  of  amplitude  modulation 
and  its  effects,  consider  the  simple  system  shown  in  Fig.  2-1.  In  this  system,  a  general  input 
signal  a{t)  is  fed  to  a  multiplier  that  mixes  it  with  another  signal,  b{t). 


a{t) 


>  d{t) 


Figure  2-1:  Mixer  Adder 
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The  sinusoidal  signals  a{t)  and  b(t).can  be  written  as  Eq.  2-la  and  2-2b, 


a(t)  =  Acos{q)J  )  (2-1) 

bit)  =  B  cos(ty^  (2.2) 

where  A  and  B  are  amplitudes,  cua  and  cob  are  angular  frequencies,  and  ^  is  a  phase  angle.  The 
multiplier  output,  which  is  denoted  as  c(0,  can  be  written  as: 


c{t)  =  ^^cos((<y^  -o)g)t  -^)  +  ^^cos((<y^  +Q}g)t  +  (p) 


(2-3) 


Following  the  multiplication  process,  a(t)  is  added  to  the  multiplier  output,  producing  the  signal. 


A.B  AB 

d(t)  =  Acos{a)^t)-\ - cos{{co^  - cOg) t  - (^)  +  — cos((<y^ -i- 

3 


(2-4) 


1 


2 


Note  that  the  output  d(t)  is  comprised  of  three  components: 

1.  A  carrier  signal,  a(t). 

2.  A  signal  at  the  difference  frequency,  coa  -  (Ob 

3.  A  signal  at  the  sum  frequency,  (Oa  +  (Ob. 


In  many  diagnostic  problems,  a  motor  current  will  take  a  form  similar  to  the  signal  dit). 
In  that  case,  coa  would  be  the  line  frequency  (i.e.  60Hz  or  12071  rad/sec)  and  cob  would  be  the 
frequency  of  some  modulating  signal.  To  detect  a  fault  in  question,  one  would  search  the 
frequency  spectrum  for  the  signals  at  the  sum  and  difference  frequencies. 

To  demonstrate  the  effects  of  amplitude  modulation,  the  system  shown  in  Fig.  2-5  was 
simulated  in  Matlab.  The  two  input  signals  are  shown  in  Fig.  2-2a  and  2-2b.  The  first  is  a  60Hz 
input  provided  as  the  carrier  signal  ait);  the  latter  is  a  lOHz  waveform  provided  as  the 
modulating  input  bit).  The  mixer  output  c(0  is  shown  in  Figure2-2c,  and  the  final  output  dit)  is 
shown  in  Fig.  2-2d.  These  plots  were  generated  using  the  constants  provided  in  Table  2-1 
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Table  2-1:  System  Inputs 


Amplitude 

Erequency 

Phase  Angle 

a{t) 

IV 

bit) 

IV 

coa 

271*60Hz 

coa 

271*10Hz 

(/> 

0  rads 

Fig.  2-3  shows  the  magnitude  of  the  frequency  spectrum  of  the  of  the  signal  d{t)  from 
Fig.  2-6.  According  to  Eq.  2-3,  this  spectrum  should  display  impulses  at  the  frequencies  Wa,  Wa- 
0B  ,  and  a)A  +  Mb.  Given  the  parameter  values  in  Table  2-1,  the  expected  frequencies  are  60Hz, 
50Hz,  and  70Hz,  respectively.  Given  that  phi  =  0  and  A=B,  Eq.  2-3  shows  that  the  amplitude  of 
the  signals  at  the  sum  and  difference  frequencies  are  exactly  one-half  the  amplitude  of  the  carrier 
signal.  These  predictions  match  the  result  shown  in  Eig.  2-3. 


Time  (sec) 


Figure  2-2a:  The  carrier  signal,  a(t).  In  this  example,  the  frequency  of  this  signal  is  60Hz. 


Figure  2-2b:  The  modulating  signal,  fe(t).  In  this  example,  the  frequency  of  this  signal  is  lOHz. 
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Time  (sec) 

Figure  2-2c:  The  mixer  output,  c{t). 


O  20  40  60  80  1 0O  1 20 


Frequency:  Hz 


Figure  2-3:  FFT  of  d{t)  from  Fig.  2-6d.  Note  that  the  amplitude  of  the  sidebands  is  approximately  one-half  the 
amplitude  of  the  60Hz  carrier  signal.  This  is  predicted  by  Eq.  2-3. 
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2.2  Using  MATLAB  to  produce  spectral  measurements 

There  are  typically  two  basic  forms  of  data  that  can  be  recorded,  stored,  and/or  displayed 
by  the  NILM  system.  Those  two  forms  are  voltage  current  pairs,  sometimes  referred  to  as  “raw” 
and  preprocessed  data.  For  this  discussion,  it  is  assumed  that  a  single  current  and  voltage 
measurement  is  taken  using  one  voltage  transducer  and  one  current  transducer  monitoring  a 
single  phase  inductive  load.  It  is  also  assumed  that  the  second  generation  NILM  software  and 
hardware  are  used,  which  is  discussed  in  chapter  5  of  this  thesis.  “Raw”  data  is  recorded  in  a 
data  channel  pair,  one  for  measured  voltage  and  one  for  current.  Power  is  calculated  with  the 
pre-processor  software,  which  essentially  multiplies  the  two  “raw”  channels  or  voltage  and 
current.  Recall  that  each  of  these  channels  for  a  given  dataset  can  be  represented  by  Eq  2-1  and 
2-2.  Preprocessed  data  takes  in  two  channels  and  produces  an  even  number  of  channels 
corresponding  to  fundamental  real  and  reactive  power  calculations  as  shown  in  figure  2-7.  Both 
“raw”  and  “preprocessed”  data  presented  in  this  thesis  was  generated  manually  using  a  simple 
MATLAB  script.  This  data  is  presented  with  an  emphasis  on  useful  spectral  quantities  and 
minimizes  the  presentation  of  non-reactive  loads  or  DC  components.  Therefore,  it  is  important 
that  the  method  of  both  “raw”  and  “preprocessed”  data  presentation  is  explained. 

The  general  Matlab  script  used  to  present  raw  current  fast  Fourier  transforms  (FFTs)  is 
presented  below  in  the  left  column  and  explained  in  the  right  column  with  examples  of  the  data 
at  significant  points  on  the  following  page.  In  this  case  the  “file”  is  a  60  Hz  signal  modulated  at 
10  Hz,  which  was  developed  using  methods  presented  later  in  this  chapter. 


%  For  "Bench  Test"  Test 
%  60Hz  signal  amplitude  1  Volt  peak  -  peak 
%  Data  taken  for  30  seconds 

f_samp  =  8000; 

Standard  sample  rate  for  NILM  is  8kHz 

load  ‘file’ 

Load  specific  files  for  test 

start  =  2000; 

Index  where  the  FFT  Starts,  use  this  value  to  move 

FFT  window  through  the  dataset.  Must  be  >  1 . 

file  E  =file(:,1); 

filej  =file(:,2); 

Define  each  data  vector  (voltage  and  Current)  (Fig  2-4) 

Ncyc  =  200; 

Number  of  cycles  used  in  the  FFT,  Influences  FFT 
amplitude. 

f_wall  =  60; 

Frequency  of  Alternating  Current. 

N_samp  =  round((N_cyc/f_wall)*f_samp); 

Normalizes  vector  length 

ind  =  (start:start+N_samp-1)'; 

Sample  pointer  initializing  at  start  variable 

fileJJ  =  file_l(ind); 

Indices  (ind)  to  be  represented  as  vectors  (Fig  2-5) 

file_l_Ldet  =  detrend(file_l_i); 

Removes  the  mean  value  or  linear  trend  from  the 
specified  dataset  vector.  (Fig  2-6) 
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hannmeister  =  hanning(N_samp); 

Creates  a  Hanning  window  to  alleviate  spectral  leakage. 

It  is  important  to  note  that  using  this  function  will  reduce 
the  amplitude  of  the  FFT  by  half. 

Nfs  =  2''16; 

Sets  the  number  of  points  in  the  FFT 

f_x_file  =  abs(fft(hannmeister.*file_l_i_det,  Nfs)); 

Contains  the  magnitude  of  the  FFT 

fx_x  =  (0:Nfs/2-1)'*(f_samp/Nfs); 

Frequency  vector  corresponding  to  the  first  half  of  the 

FFT 

plot  (fx_x,f_x_file(1  :length(fx_x)),'b'); 

Plotting  the  resulting  FFT  (Fig  2-7) 

Channel  Data 


Samples  (@8000  eye  /  see) 

Figure  2-4:  Sample  Dataset 


Dynamic  range  of  the  analog  to  digital  device  is  4096  amplitude  divisions.  Thus,  data 
points  above  2048  are  positive  and  points  below  2048  are  negative  with  2048  representing  zero. 
Sample  dataset  presented  in  figure  2-5  oseillates  around  an  amplitude  or  “y”  value  of  2048. 
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Segment  of  Data  to  be  Processed 


Figure  2-5:  Part  of  Signal  used  in  FFT 


Channel  Detrend 


Figure  2-6:  Detrending  selected  dataset 
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FFT  of  60Flz  Modulated  by  10Hz 
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Figure  2-7:  Resulting  FFT  in  Matlab 


It  can  be  shown  that  Fig.  2-3  and  2-7  are  the  same  with  the  exception  of  the  50  % 
reduction  in  amplitude  caused  by  using  the  Hanning  window  [19]. 


2.3  The  NILM  Preprocessor 

The  NILM  collects  and  processes  two  input  data  streams  -  line  voltage  and  aggregate 
current.  Using  these  two  quantities,  software  on  the  NILM  can  perform  a  number  of  important 
tasks,  including: 

•  Real-power  computation 

•  Reactive-power  computation 

•  Current-harmonic  computation 

•  Load  identification 

•  Diagnostic  analysis 

•  Power-quality  analysis 

The  first  three  tasks  listed  above  are  performed  by  the  NILM’s  preprocessor.  This 
software  package,  which  is  commonly  referred  to  as  “prep,”  performs  a  spectral  decomposition 
of  the  current  by  computing  estimates  of  the  quantities  defined  in  Eq.  1-1  and  1-2.  This  section 
proves  a  qualitative  description  of  the  tasks  performed  by  prep,  and  it  also  describes  how  the  user 
should  interpret  various  preprocessed  outputs  in  typical  monitoring  situations. 
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2.3.1  Description  of  Operation 


Figure  2-8  presents  a  block  diagram  of  the  preprocessor.  As  shown,  prep  takes  two 
inputs  -  line  voltage  and  measured  current  -  and  it  produces  eight  outputs.  The  exact  meaning  of 
each  output  depends  on  the  specific  details  of  a  particular  installation.  Later  sections  will 
describe  two  typical  installations. 


System  Voltage 
System  Current 


ai,bi 


as.ba 


a5,b5 


a?,  by 


Figure  2-8:  Preprocessor  Block  Diagram 


Assuming  minimal  distortion,  the  line  voltage  can  be  approximated  as  either  a  pure 
cosine  or  as  a  pure  sine.  In  this  example,  it  will  be  assumed  that  the  line  voltage  is  a  cosine,  Eq. 


2-4. 


v(t)  =V  cos{cot) 


(2-5) 


In  general,  the  current  drawn  from  such  a  source  will  contain  a  number  of  harmonics,  i.e. 


i(t)  -  /j  cos  {cot  +  )-i-  /j  cos  {icot  -I-  (Z>2 )+  /3  cos  {fcot  -I-  ^3  )-i-  ••• 

=  ^  cos  {kcot  +  (p^ ) 

k 


(2-6) 


Equation  2-2  highlights  a  key  fact  -  the  current  drawn  at  any  one  harmonic  frequency  is 
characterized  using  two  quantities,  namely  an  amplitude  h  and  a  phase  shift  ^k.  The  phase  shift 
represents  any  offset  from  the  line  voltage,  and  the  amount  of  phase  shift  is  strongly  dependent 
on  the  type  of  load.  Purely  resistive  devices  such  as  heating  elements  draw  currents  that  are 
completely  in  phase  with  the  voltage,  i.e.  (pi  =  0.  By  comparison,  a  purely  inductive  or 
capacitive  load  will  draw  a  current  that  is  90  degrees  out-of-phase  with  the  line  voltage,  i.e.  (p  i  = 
+90°.  These  loads  are  often  referred  to  as  purely  reactive.  Many  common  devices  such  as  motors 
consist  of  both  resistive  and  reactive  elements.  In  those  cases,  the  magnitude  of  the  phase  shift 
will  lie  between  0  degrees  and  90  degrees.  The  sign  of  the  phase-shift  depends  on  whether  the 
load  is  either  inductive  or  capacitive  in  nature. 
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Equation  2-5  suggests  that  loads  draw  currents  that  are  rich  in  harmonics.  This  general 
statement  reflects  the  fact  that  loads  can  be  non-linear.  Examples  of  non-linear  loads  include 
devices  such  as  computers,  motor  controllers,  and  fluorescent  lamps.  Typically,  these  non-linear 
devices  contain  a  power  supply  with  a  rectifier.  As  an  example,  Eig.  2-9  presents  the  current 
drawn  by  a  small  personal  computer.  As  shown,  the  current  is  neither  a  pure  sine  nor  a  pure 
cosine.  As  predicted  by  system  theory,  the  current  drawn  by  the  device  can  be  represented  as  a 
Eourier  series  of  the  form  presented  in  Eq.  2-5  [20].  By  comparison,  the  currents  drawn  by 
devices  such  as  motors  and  heating  elements  draw  minimal  harmonic  content.  Thus,  these 
currents  can  be  approximated  as: 


i(t)  =  4  cos(a)t  -h  0i) 


(2-7) 


Figure  2-9:  Typical  Personal  Computer  Current  Measurement  [22] 


Figure  2-10:  Typical  Induction  Motor  Start-up  Current  Measurement  [22] 
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Figure  2-10  shows  the  current  drawn  by  a  typical  induction  motor.  Note  that  the  current 
can  be  approximated  as  a  sinusoid  with  a  time-varying  amplitude. 

It  should  be  noted  that  the  current  presented  in  Eq.  2-5  can  be  re-arranged  through  a 
series  of  trigonometric  manipulations.  In  effect,  the  current  at  each  harmonic  can  be  written  as 
the  sum  of  a  sine  wave  and  a  cosine  wave,  i.e. 

i(t)  =  /j  cos(^j  )cos(<J>r) -  /j  sin((Zij  )sin(<af) -I- cos((Zi2  )cos(2<af) - sin(^2  )sin(2<j>r)  H — 

The  cosine  at  each  frequency  represents  the  resistive  component  of  the  current  and  the 
sine  represents  the  reactive  component.  Thus,  purely  resistive  devices  draw  currents  containing 
only  cosine  terms  and  purely  reactive  devices  draw  currents  containing  only  sine  terms.  Loads 
with  both  resistive  and  reactive  elements  draw  a  combination  of  both.  The  representation 
presented  in  Eq.  2-7  will  be  useful  in  the  discussion  to  follow. 

The  primary  function  of  the  preprocessor  is  to  estimate  the  real  and  reactive  power  drawn 
from  the  ideal  single-phase  source  presented  in  Eq.  2-1.  By  definition,  the  real  power  is 

Ej  =  —  [  V  cos{ax)i{t)dt 

T  (2-9) 

and  the  reactive  power  is 


2i  ^^\Vsm{Qx)i{t)dt 
I 


(2-10) 


If  i{t)  is  of  the  form  presented  in  Eq.  2-6,  then  these  results  simplify  to 


)cos"(6;f>/t  = 

J  jj  2 


(2-11) 


and 


sin((Z>,  )un^{Qx)dt  =  - ^ 

T  JT  2 


(2-12) 


Similarly,  prep  also  estimates  the  in-phase  and  quadrature  components  of  the  current  at 
several  odd  harmonics.  At  the  third  harmonic,  for  instance,  it  estimates  the  following  two 
quantities: 


E3  =  ^|^cos(3tyr)?(()c/(  =  ^3  (^3 


(2-13) 


25 


(2-14) 


03  =  ^|^sin(3<2lf)/(r)(ir  = 

Similar  calculations  are  currently  performed  at  both  the  5*  and  7*  harmonie,  although 
others  eould  be  eonsidered  if  desired. 


The  eomplete  details  of  the  power  ealeulations  performed  by  the  preproeessor  are 
relatively  intensive,  and  they  are  beyond  seope  of  this  thesis  [3].  For  purposes  of  illustration,  it 
oan  be  assumed  that  prep  generates  a  eolleetion  of  sinusoids  internally.  As  presently 
implemented  [20],  the  following  signals  are  ereated: 


=  cos(<si),  v^j  =  sin(<si),  v^3  =  cos(3tyt),  v^3  =  sin(3<si), 
v^5  =  eos(5<af),  =sin(5<Mf),  =eos(7ty!f),  v^-j  -smijox) 


These  phase-loeked  waveforms  are  used  to  estimate  the  power-related  quantities 
discussed  above.  For  example,  the  preprocessor  uses  Vai  to  estimate  the  following  integral: 


(2-16) 


This  result  can  be  related  to  P/  if  sealed  by  the  amplitude  of  the  measured  line  voltage. 
Similar  results  hold  for  the  other  eoeffieients. 


Because  the  amplitude  and  phase  of  the  eurrent  ehange  over  time,  the  output  of  the 
preproeessor  must  ehanges  as  well.  The  preproeessor  eomputes  the  eoeffieients  Uk  and  bk  over  a 
moving  window.  This  window,  whieh  is  synchronized  to  the  measured  line  voltage,  is  advaneed 
every  half-cyele.  Thus,  new  eoeffieients  are  produeed  at  twice  the  line  frequeney.  For  a  60Hz 
service,  this  implies  that  power  estimates  are  produeed  at  120Hz. 


After  being  ealeulated  by  the  preproeessor,  the  eoeffieients  ai  and  bi  can  be  sealed  by  the 
amplitude  of  the  line  voltage  in  order  to  produee  estimates  of  Pi  and  Qi.  These  values  will  vary 
slowly  over  time  and  can  be  used  for  identifieation  purposes.  Exact  details  of  the  eomputation 
process  are  provided  in  [1]. 
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2.4  Bench  Testing 


The  first  step  in  the  proeess  in  eonfirming  the  NILM  is  eapable  of  deteeting  faults  that 
other  vibration  analysis  and  eurrent  measurement  analysis  (CMA)  teehniques  ean  is  to  present 
the  NILM  with  a  realistic  known  input  and  compare  the  output  with  the  predicted  output.  Figure 
2-11  illustrates  the  conceptual  approach  to  developing  this  test. 
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Figure  2-11:  Block  diagram  of  Fundamental  Test 


2.5  Test  Signal  Generation 

For  this  test  setup,  the  Tektronix  AWG2005  arbitrary  waveform  generator  was  used  to 
generate  two  waveforms  representing  current  and  voltage.  CMA  reveals  that  a  multitude  of 
mechanical  modes  in  induction  machines  can  be  monitored.  Nearly  every  electromechanical 
system  has  at  least  one  significant  mode.  Items  such  as  belts,  fans,  impellors,  and  pumps  exhibit 
modes  that  can  be  monitored  in  a  demodulated  current  spectrum  at  a  frequency  representative  of 
the  physical  system.  To  create  a  test  that  emulates  the  simplest  scenario  for  an  electromechanical 
system  a  60Hz  waveform,  representing  supply  voltage  and  an  Amplitude  Modulated  (AM) 
waveform  representing  current  with  a  single  mechanical  mode  was  created.  A  sample  rate  of  60 
kHz  was  used  to  generate  each  of  these  waveforms.  In  order  to  generate  the  AM  waveform,  the 
AWG2005  requires  a  separate  input  to  modulate  the  base  frequency  (60  Hz  in  this  case).  A 
Hewlet  Packard  8656B  signal  generator  was  used  as  the  modulating  input  to  the  AWG2005. 
AWG2005.  Two  channels  of  the  AWG2005  and  an  oscilloscope  were  used  to  ensure  proper 
amplitude,  frequency,  and  phase  was  entering  the  NILM  box.  Figure  2-12  on  the  following  page 
presents  the  test  setup  described  above  in  LEES. 
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Figure  2-12:  "Bench  Test"  setup 


Channel  1  of  the  Tek  AWG2005  is  the  only  channel  that  can  be  modulated.  Therefore 
the  simulated  current  originates  here  but  enters  channel  2  of  the  Oscilloscope.  This  was  done 
because  the  voltage  will  be  the  reference  for  current  phase  shifting.  Therefore,  this  set-up  will 
make  phase  adjustments  more  straightforward  in  subsequent  test. 

In  initial  NILM  installations,  sensors  are  connected  directly  up  stream  of  the  load  and 
tailored  to  that  specific  load.  The  voltage  sensor  selected  in  most  application  is  the  LV  25 -P 
manufactured  by  LEM.  The  current  sensor  selection  requires  a  closer  look  at  both  the  transient 
current  signatures  as  well  as  typical  steady  state  current  ranges  within  expected  load  variations. 
Usually,  a  current  sensor  is  selected  based  on  the  maximum  expected  current  spike  of  a  particular 
load.  The  sensor  is  then  scaled  as  a  function  of  NILM  gain,  selected  current  transducer  (CT) 
conversion  ratio  and  measuring  resistor(s).  In  general,  current  transducers  are  ideal  current 
sources.  The  magnetic  flux  created  by  the  primary  current  is  balanced  by  a  complementary  flux 
produced  by  driving  a  current  through  the  secondary  windings.  A  hall  device  and  associated 
electronic  circuit  are  used  to  generate  the  secondary  (compensating)  current  that  is  an  exact 
representation  of  the  primary  current.  Thus,  the  single  conductor  output  of  the  CT  represents  an 
ideal  current  source  representation  of  the  current  passed  through  the  CT  [21]. 

The  simplest  NILM  installation  uses  two  input  channels,  one  channel  to  monitor  load 
voltage  and  one  to  monitor  load  current.  For  this  first  test,  we  used  two  voltage  sources,  the  first 
representing  the  voltage  CT  input  and  the  second  voltage  source  was  passed  across  a  lOOQ 
resistor  to  develop  a  current.  Figure  2-14  shows  a  screen  capture  of  the  oscilloscope  with  the 
first  test  signal  driven  to  the  NILM  box.  It  presents  both  the  “voltage  waveform”  (in  blue  green) 
and  the  “current  waveform”  (in  yellow).  Ideally,  this  represents  an  electromechanical  system 
that  has  one  electro  mechanical  mode.  In  this  case  a  current  modulation  frequency  of  10  Hz  was 
used.  The  amplitude  of  the  simulated  voltage  is  1  volt  peak-to-peak  and  the  modulated  voltage 
or  “current”  is  1  volt  peak-to-peak.  This  was  done  to  provide  enough  separation  that  scaling 
factors  in  previous  studies  like  the  one  done  by  Patrick  Bennett  [20]. 
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Figure  2-13:  Oscilloscope  screen  capture  of  both  the  60  Hz  “voltage”  waveform  and  60  Hz  “current” 

waveform  modnlated  at  10  Hz  in  phase 


Fignre  2-14:  Oscilloscope  captnre  of  both  60  Hz  "voltage”  and  60  Hz  current"  waveform  modnlated  at  10  Hz 

and  leading  voltage  by  45  degrees 


The  modulated  signal  does  not  appear  to  quite  make  1  V  pp  in  either  Fig.  2-13  or  2-14. 
Flowever,  these  amplitudes  were  direet  data  entries  into  the  AWG2005  user  interfaee;  therefore, 
no  effort  was  made  to  eompensate. 


The  signals  above  were  fed  direetly  into  the  terminal  bloek  seen  in  Fig.  2-15  on  the 
following  page.  The  waveforms  were  relatively  small  eompared  to  standard  shipboard 
applieations  to  demonstrate  the  NILM’s  ability  to  handle  small  signals  and  observe  system  noise. 
Figure  2-15  shows  the  signal  entering  the  NILM  input  terminal  bloek. 
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Figure  2-15:  Test  Signal  Generator  output  to  NILM  connection 


2.6  NILM  Configuration  and  Measurement 

The  NILM  computer  was  configured  similar  to  standard  marine  installation  with  a  couple 
of  exceptions.  In  order  to  look  at  the  Analog  to  Digital  A/D  gain  function,  the  gain  code  for  the 
voltage  channel  was  set  to  1,  meaning  it  would  appear  to  be  doubled.  The  gain  code  for  the 
current  channel  was  set  to  zero  for  a  1-to-l  relationship.  The  NILM  voltage  and  current 
channels  were  sampled  at  a  rate  of  8  kHz  for  30  seconds. 


2.6.1  Representing  Current 

In  order  to  generate  a  clean  realistic  signal,  a  waveform  generator  (TEK  AWG2005)  was 
used  which  imposes  some  minor  modeling  limitations.  Most  importantly,  the  output  of  the 
function  generator  more  closely  resembles  an  ideal  voltage  source  whereas  the  output  of  the 
typical  LEM  current  transducer  output  looks  more  like  an  ideal  current  source.  The  output  of  the 
waveform  generator  was  driven  across  a  resistance  inducing  a  representative  current.  However, 
to  stay  true  to  the  output  of  the  test  setup,  both  channel  1  and  2,  representing  voltage  and 
“current”  respectively,  will  both  be  referred  to  as  voltages.  Channel  2  instantaneous  voltage  can 
be  divided  by  the  terminal  board  resistance  to  determine  NILM  current.  This  will  also  be 
important  in  representing  power. 

Both  the  voltage  and  current  signals  were  driven  into  a  lOOQ  resistor  (across  the  terminal 
board  pins  33  to  67  and  34  to  68  respectively).  Comparing  figures  2-7  and  2-3,  one  can  see  that 
the  NILM  has  effectively  captured  the  exact  input  waveform  created  by  the  test  setup.  Referring 
back  to  figure  2-13,  the  agreement  between  the  oscilloscope  depiction  of  the  input  and  figure  2- 
16  on  the  following  page  confirm  NILM  operation  up  to  the  preprocessor. 
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2-16:  Channel  Output  from  the  NILM 


In  general,  the  amplitude  of  the  voltage  representing  a  current  for  such  a  mixed  signal  can 
be  calculated  by  Eq.  2-16  where  A  is  the  amplitude  of  the  carrier,  M  is  the  amplitude  of  the 
modulation,  Mc  is  frequency  of  the  carrier  (in  radians),  Wm  is  the  frequency  of  the  modulation  (in 
radians),  t  is  time  in  seconds  and  0  is  the  phase  shift  between  the  carrier  and  modulated  signals 
(in  radians). 

To  confirm  that  there  are  no  other  “hidden”  signals  in  the  output  waveform  a  fast  Fourier 
transform  (EFT)  was  performed  on  the  raw  data  using  MATLAB.  The  EFT  should  show 
amplitude  consistent  with  the  time  domain  plot.  In  general,  the  amplitude  of  the  voltage 
representing  a  current  for  such  a  mixed  signal  can  be  calculated  by  Eq.  2-1  where  A  is  the 
amplitude  of  the  carrier,  M  is  the  amplitude  of  the  modulation,  Wc  is  frequency  of  the  carrier  (in 
radians),  Wm  is  the  frequency  of  the  modulation  (in  radians),  t  is  time  in  seconds  and  cp  is  the 
phase  shift  between  the  carrier  and  modulated  signals  (in  radians). 


Vout  =  4  cos(mct)  +  ^[(cos((mc  -  +  0)  +  cos  (((rUc  +  +  o)] 


(2-17) 
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2.6.2  Power 


The  NILM  uses  measurements  of  voltage  and  eurrent  to  produee  estimates  of  real  and 
reaetive  power  eonsumption  (heneeforth  referred  to  as  P  and  Q,  respeetively).  For  a  single  phase 
system,  the  NILM  is  supplied  with  the  voltage  from  line  to  neutral  and  the  eurrent  flowing  to  any 
loads  downstream  of  the  monitoring  point.  Real  power  (P)  is  produeed  by  eurrent  whieh  is  in- 
phase  with  voltage  and  reactive  power  (Q)  is  produced  by  current  which  is  90  degrees  out  of 
phase  with  voltage.  It  also  estimates  the  higher  harmonic  content  of  current  to  produce  the  first 
three  odd  harmonics  of  real  and  reactive  power  (referred  to  as  Pk  and  Qk  for  the  kth  harmonic  of 
P  and  Q,  respectively).  If  the  same  sinusoidal  voltage  is  applied  to  a  resistor,  capacitor,  and 
inductor,  the  NILM  would  produce  positive  P  with  zero  Q  for  the  resistor,  positive  Q  with  zero  P 
for  the  inductor,  and  negative  Q  with  zero  P  for  the  capacitor.  The  convention  of  using  a  power 
angle  measured  in  the  clockwise  direction  from  the  real  axis  is  employed,  so  that  the  power  angle 
for  an  inductor  is  -i-90  degrees  while  the  power  angle  for  a  capacitor  is  -90  degrees  [17].  For  a 
more  in— depth  description  of  the  NILM  calculations  used  to  determine  P  and  Q,  see  [20]. 

The  NILM  was  given  one  signal  pair  from  the  waveform  generator.  Table  2-1  below  lists 
the  variables  for  each  test.  All  test  used  60  Hz  as  the  principle  frequency.  Channel  2  was  chosen 
to  lead  channel  1  simply  because  current  generally  leads  voltage  in  predominantly  inductive 
circuits.  These  tests  were  run  and  preprocessed  data  recorded  for  each. 


Table  2-2:  Phaser  Test 


Test 

Modulation 
(in  Hz) 

Phase  shift 
(in  degrees) 

1 

10 

0 

2 

10 

45 

3 

20 

0 

4 

20 

45 

5 

30 

0 

6 

30 

45 

In  order  to  check  the  preprocessor  or  prep  software,  the  initial  input  waveform  depicted 
in  Fig.  2-7  was  used.  Data  from  test  number  1  (60  Hz  voltage  channel  and  the  10  Hz  modulated 
current  channel)  was  used  and  NILM  outputs  were  plotted  in  the  Fig.  2-18  and  2-19  below. 

The  fact  that  this  test  used  an  ideal  voltage  source  to  develop  a  current  across  a  resistance 
becomes  very  important  in  looking  at  representative  Power.  For  this  reason,  we  will  scale  power 
to  an  average  value  of  1.  Because  the  amplitudes  can  always  be  scaled,  the  goal  of  the  specific 
test  was  not  to  demonstrate  amplitude  accuracy  but  to  show  that  frequency  and  phase  response 
were  accurate.  However,  the  amplitude  should  reflect  a  value  that  can  be  linearly  scaled  to 
ensure  accuracy  over  the  target  spectrum. 
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P  Prep 


2-17:  Preprocessor  Output  of  Real  Power 


Q  Prep 


2-18:  Preprocessor  Output  of  Reactive  Power 


Figures  2-17  and  2-18  reflect  expected  outputs.  The  real  power  fluctuates  around  1  Watt 
while  reactive  power  fluctuates  around  zero.  Ideally  reactive  power  would  have  no  fluctuations 
but  some  can  be  seen  in  above.  However,  the  fluctuations  are  less  than  0.04  VARs  which  is  less 
than  5%  of  the  average  wattage.  These  minor  fluctuations  could  be  caused  by  preprocessor 
filtering  techniques,  internal  noise,  external  noise  or  a  little  of  all  three.  Keeping  in  mind  that  the 
input  signal  amplitudes  are  extremely  small  compared  to  typical  values,  these  fluctuations  are 
negligible. 

In  order  to  establish  a  reference,  the  expected  of  a  strictly  resistive  test  was  used  to 
normalize  the  remaining  section  FFT  plots  to  0  dB  =  250mW.  This  was  done  by  creating  a 
synthetic  voltage  and  current  signal  using  Matlab  with  a  known  amplitude  and  running  it  through 
the  preprocessor.  This  established  a  reference  for  later  test  that  adjust  the  phase  angle  between 
the  voltage  and  current.  One  note  with  regard  to  the  FFTs  of  power  in  this  section,  the  same 
Matlab  code  presented  in  the  section  2-1  was  used  to  create  these  FFTs.  The  cursor  feature  in 
the  plot  function  was  used  to  present  the  amplitude  of  the  key  spectral  quantities. 
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Real  Power  (Test  1) 


0  10  20  30  40  50  60 

Frequency:  Hz 

2-19:  FFT  of  Preprocessed  Real  Power 


CD 

■D 


2-20:  FFT  of  Preprocessed  Reactive  Power 

Note  the  plots  in  Fig.  2-19  and  2-20  are  eonsistent  with  corresponding  time  domain  plots 
on  the  previous  page.  There  is  also  another  spectral  component  that  exists  at  70Hz.  However, 
because  of  the  Nyquist  rate,  the  sampling  rate  of  prep  (120  Hz)  limits  the  frequency  of  the  results 
to  60  Hz  so  higher  spectral  quantities  due  exist  but  are  not  displayed  without  changing  the 
preprocessor. 
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Real  Power  (Test  1 ) 


E  -60 


-100 


-120 


-140 


10Hz  with  detrend 
1 0Hz  without  detrend 


6  8  10  12  14  16  18  20 

Frequency:  Hz 


Figure  2-21:  Figure  2-21  with  plot  of  lOHz  signal  without  “detrend”  function  showing  DC  component 


A  close  look  at  the  aetual  DC  component  reveals  an  amplitude  of  -Fl2db.  This  could  be 
caused  by  test  equipment  coupling  or  preprocessor  filtering  techniques. 
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2.6.3  Applying  a  45  degree  leading  phase  shift  to  “Simulated  Current” 

The  waveform  generator  was  adjusted  to  perform  the  second  test.  This  was  done  by 
advancing  the  phase  of  the  signal  being  driven  into  channel  two  of  the  NILM  by  45  degrees. 
Figure  2-22  presents  the  scope  display  capture  of  channel  one  (blue  green)  and  channel  two 
(yellow)  fed  into  the  NILM  terminal  board.  The  only  adjustment  made  to  the  test  setup  in  test 
one  was  the  adjustment  of  the  phase  angle  between  channel  one  and  two  (channel  two  now  leads 
channel  one  by  45  degrees). 


2-22:  Oscilloscope  capture  of  60Hz  "voltage"  and  60  Hz  modulated  at  10  Hz  "simulated  current" 


10  Hz  Modulation  Leading  Voltage  by  45  Degrees 


As  the  modulated  voltage  representing  current  was  shifted,  the  amplitude  appeared  to 
decrease.  The  FFT  of  the  prepped  data  will  present  a  better  picture  of  what  is  actually  happening 
here. 
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Next,  the  phase  shifted  “simulated  eurrent”  seen  in  Fig.  2-16  was  fed  into  the 
preproeessor  software.  The  preproeessor  output  in  Fig.  2-24  and  2-25  (real  and  reactive  power 
respectively)  appear  expectedly  different. 


P  Prep 


2-24:  Preprocessor  Output  of  Real  Power  (Simulated  Current  Leads  Voltage  45  deg) 


Q  Prep 


2-25:  Preprocessor  Output  of  Reactive  Power  (Simulated  Current  Leads  Voltage  45  deg) 

Observing  Fig.  2-24  and  2-25,  it  appears  that  the  preprocessor  software  is  reflecting  the 
fact  that  both  real  and  reactive  power  shows  up  on  a  DC  level  that  is  equal  to  lW(2),  which  is 
what  would  be  expected  with  a  45  degree  shift  between  voltage  and  current.  A  look  at  the  FFTs 
of  real  and  reactive  power  should  confirm  the  time  domain  data. 
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Figure  2-26:  Preprocessor  Outputs  (Channels  1  and  2) 


Figures  2-26  is  presented  simply  to  show  an  amplitude  plot  of  the  entire  preprocessor 
output.  In  order  to  get  more  value,  Fig  2-27  focuses  in  on  areas  around  10  Hz. 


Real  Power  Prep  Result 


Figure  2-27:  Closer  look  at  Preprocessor  Real  Power 


The  unshifted  lOHz  modulation  spectral  component  shows  up  at  lOHz  and  Minus  3  dB. 
The  spectral  component  shifted  45  degrees  shows  a  1.5dB  reduction  in  amplitude  which  again  is 
expected.  A  45  degree  shift  can  be  shown  to  reduce  the  DC  level  of  the  spectrum  to  1/^2. 
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Reactive  Power  Prep  Result 


Frequency:  Hz 

Figure  2-28:  Closer  look  at  Preprocessor  Reactive  Power 

There  are  two  interesting  things  about  Fig.  2-28.  First  is  that  the  reaetive  power  in  the 
phase  shifted  lOHz  modulated  ease  does  in  fact  show  a  minus  4.5dB  amplitude.  And  second,  the 
inherent  noise  of  the  system  can  be  seen.  We  know  there  should  be  little  to  no  reactive 
component  in  the  unshifted  case.  Fig.  2-28  indicates  that  a  component  with  an  amplitude  of 
approximately  -23db,  which  translates  to  less  than  a  femtowatt.  The  remaining  test  provided 
similar  results  on  all  accounts. 


2.6.4  Conclusion 

Electromechanical  systems  are  inherently  designed  with  mechanical  resonating  modes 
that  can  be  detected  in  electrical  current  signatures.  This  is  one  of  the  fundamental  principles  the 
NILM  system  exploits.  The  NILM  records  and  processes  voltage  and  current  relationships  and 
presents  robust  data  sets,  including  but  not  limited  to  real  and  reactive  power  relationships,  that 
can  be  used  to  make  technical  decisions.  The  tests  above  sought  to  use  a  two  channel  waveform 
generator  to  model  the  simplest  electromechanical  relationship,  an  electromechanical  system 
with  a  single  resonate  mode  with  the  goal  of  observing  the  NILMs  accuracy  in  data  recording, 
multi-channel  phase  and  amplitude  relationships,  preprocessor  operation  and  the  level  of  system 
noise.  This  simple  test  confirmed  the  standard  NILM  setup  operates  as  designed  in  each  of  the 
observed  cases. 
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3  AC  Motor  Analysis 

3.1  Introduction 


The  NILM  is  able  to  monitor  and  report  on  a  number  of  key  motor  faults.  The  NILM 
performs  these  diagnosties  using  both  the  eurrent  spectrum  and  the  spectral  envelopes  discussed 
previously.  The  use  of  the  NILM  in  this  capacity  is  extremely  attractive  because  the  NILM  itself 
is  installed  in  a  central  location  and  does  not  require  special  motor  access.  Additionally,  the 
NILM  offers  the  opportunity  to  consider  monitoring  multiple  motors  simultaneously. 

The  next  two  chapters  provide  an  initial  demonstration  of  the  NILM’s  capabilities  as  a 
motor-monitoring  device.  This  chapter  begins  by  examining  the  use  of  the  NILM  to  monitor  a 
single  motor;  Chapter  4  considers  the  possibility  of  monitoring  multiple  motors  by  applying  the 
process  used  in  this  chapter  but  using  an  aggregate  current  signal. 


3 . 2  General  Procedure 

Using  information  contained  in  [24],  the  most  relevant  motor-related  problems  were 
identified  and  bulleted  on  the  following  page.  The  two  most  common  problems,  motor 
imbalance  and  bearing  faults,  drove  the  development  of  specific  tests  in  laboratory.  These  test 
involved  imposing  faults  representative  of  rotor  imbalances  and  bearing  faults  seen  in  the  field. 
The  ultimate  vision  of  this  application  is  that  an  aggregate  sensor  will  be  used  to  monitor 
multiple  machines  in  an  automated  fashion.  This  vision  shapes  the  way  the  NILM  is  configured 
to  monitor  and  the  procedure  used. 

Ultimately,  the  goal  is  to  develop  a  method  to  autonomously  monitor  and  report  on 
multiple  motors  with  a  single  aggregate  sensor.  With  this  in  mind,  it  is  clear  that  a  method  must 
be  established  to  autonomously  monitor  and  report  on  a  single  motor  with  a  single  sensor.  A 
process  is  presented  that  involves  defining  physical  machine  characteristics,  monitoring  the 
current  spectrum  to  confirm  a  machine  with  defined  characteristics  can  be  identified,  “locking 
on”  to  spectral  current  components  associated  with  that  machine  to  track  speed  and  /  or  slip,  and 
using  that  speed  /  slip  with  the  corresponding  spectral  components  associated  with  imbalance  or 
bearing  faults  to  determine  machine  health. 


3.3  Common  Motor-Related  Problems 

The  diagnostics  discussed  in  this  chapter  involve  either  a  spectral  decomposition  of  the 
current  or  its  spectral  envelopes  (i.e.  Pi  and  Qi  from  Chapter  2).  In  order  to  develop  a  NILM 
diagnostic  system  to  determine  the  health  of  motor  driven  systems,  some  specific  initial  target 
problems  should  be  identified.  As  a  starting  point,  key  faults  that  were  known  to  affect  the 
current  spectrum  [15]  were  investigated;  they  are  listed  on  the  following  page. 
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•  An  unbalanced  rotating  element 

•  Bearing  faults 

•  Misalignment  between  motor  and  load 

•  Bent  shafts 

•  Shaft  radial  thermal  expansion 

•  Cavitation 

•  Water  hammers 

•  External  vibration 

The  first  two  bullets  above  represent  the  most  commonly  documented  root  causes  of  motor 
failure  [24],  Initially,  a  simple  single  phase  axial  fan  motor  were  used  to  establish  the  NILM’s 
ability  to  detect  evidence  of  an  imbalanee  and  bearing  faults. 

3.4  Developing  NILM-based  Diagnostics  for  Induction  Machines 

To  detect  faults  such  as  the  ones  discussed  previously,  a  NILM  monitoring  aggregate 
current  flow  must  follow  a  systematic  procedure  consisting  of  both  load  identification  and 
subsequent  diagnosis.  This  section  describes  the  process  that  has  been  developed  during  the 
course  of  this  project.  The  approach  has  been  used  to  monitor  both  single  and  multiple  motor 
loads.  Figure  3-1  shows  a  block  diagram  of  the  proposed  monitoring  process.  Digital  stator 
current  data  is  passed  to  a  signal  processing  routine,  which  selects  and  captures  a  specific 
window  of  time  domain  data.  The  fast  Fourier  transform  is  generated  and  the  resulting  spectral 
current  data  is  first  used  to  confirm  the  presence  of  induction  machine  loads.  For  this  specific 
example,  motor  slot  harmonics  are  identified  using  key  motor  parameters  entered  into  a  graphical 
user  interface  (GUI).  This  GUI  should  be  similar  to  the  one  designed  in  [23]  but  would  present 
current  diagnostic  information.  The  GUI  would  have  a  data  entry  function  where  the  key  motor 
parameters  for  the  machines  under  test  would  be  entered.  This  information  would  be  fed  into  the 
Information  and  Computation  function,  which  stores  the  key  parameters  and  uses  a  slot  harmonic 
calculation  to  determine  slip.  It  then  packages  this  information  in  a  vector  and  passes  it  to  the 
diagnostic  subroutines  according  to  preferences  entered  into  the  GUI. 


Stator 

Current 

Data 


Figure  3-1:  General  Single  Motor  Diagnostic  Process  Block  Diagram 
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The  diagnostic  modules  are  arranged  in  a  hierarchy  associated  with  the  needs  of  the  end 
user.  The  output  to  each  loaded  diagnostic  is  controlled  by  a  graphical  user  interface  based  on 
system  requirements.  There  are  some  steps  in  this  process  that  need  to  be  discussed.  Those 
steps  are  presented  in  the  following  sub-section. 


3.4.1  Identify  key  motor  parameters 

The  first  step,  which  must  be  undertaken  before  monitoring,  is  to  determine  certain  key 
motor  parameters  such  as  the  line  frequency,  the  number  of  pole  pairs,  bearing  geometry, 
machine  loading,  and  environment.  Most  of  the  machine  parameters  needed  to  make  initial  slip 
calculations  can  be  found  on  the  machine  name  plate.  Given  that  ultimately  an  aggregate  sensor 
will  be  used  to  monitor  multiple  loads,  a  single  system  supply  frequency  can  be  assumed.  In  this 
case,  the  current  spectrum  FFT  can  be  used  to  determine  the  system  supply  frequency  by  looking 
at  the  spectral  component  with  the  highest  amplitude  and  confirming  this  using  supply 
harmonics.  Additional  information  required  can  be  obtained  from  the  manufacturer,  such  as  the 
number  of  rotor  slots,  bearing  geometry,  etc.  The  complete  list  of  key  parameters  is  based  on  the 
type(s)  of  diagnostics  desired. 


3.4.2  Estimate  speed 

The  second  step  in  the  monitoring  process  is  to  identify  motor  speed.  This  requires  a 
combination  of  NILM’s  powerful  load-classification  routines  and  the  detection  of  certain  key 
spectral  quantities.  The  NILM’s  load  classifier,  which  is  described  in  [23],  can  identify  specific 
loads  operating.  Once  such  a  decision  has  been  made,  the  NILM  can  use  the  parameters 
determined  previously  to  estimate  the  speed  of  that  motor. 

A  number  of  sensorless  speed  estimation  procedures  have  been  identified  over  the  years 
[24].  The  estimation  procedure  used  in  this  thesis  relies  on  the  successful  identification  of 
several  rotor-slot  harmonics.  Further  information  about  these  signals  and  the  speed-estimation 
procedure  are  presented  in  Section  3.4.3.  If  identical  motors  are  operating,  the  rotor-slot 
harmonics  could  be  used  to  assist  in  the  identification  process.  More  details  are  presented  in 
chapter  4. 


3.4.3  Verification  of  Load  Identification  and  Speed  Estimate 

Rotor-slot  harmonics  are  not  the  only  speed-related  signals  present  in  the  current 
spectrum.  One  other  speed  indicator  that  can  be  used  to  verify  the  speed  estimate  is  eccentricity- 
related  components  that  appear  as  modulations  of  the  60Hz  line  current.  These  signals  can  also 
be  used  to  detect  certain  critical  faults  such  as  load  imbalances.  While  these  spectral  components 
are  not  as  accurate  as  slot  harmonics,  they  are  valuable  in  the  role  of  confirming  the  primary 
speed  estimation  method. 


42 


3.4.4  Health  Monitoring  and  Fault  Detection 


Once  motors  have  been  successfully  identified,  it  is  possible  to  start  tracking  certain  key 
spectral  quantities  that  are  indicative  of  failure.  In  this  thesis,  two  different  types  of  faults  are 
considered  -  load  imbalances  and  bearing  wear.  The  first  of  these  can  be  monitored  using  the 
eccentricity-related  signals  described  previously;  the  latter  can  be  monitored  by  tracking  certain 
characteristic  bearing  frequencies  that  are  present  in  the  current  spectrum.  Further  details  are 
presented  in  Section  2. 6. 2. 3. 


3.5  Experimental  Verification 

The  steps  in  the  process  described  in  Section  3.4  were  experimentally  verified  in  the 
laboratory.  Initial  experiments  were  performed  using  a  single-phase  axial  fan  motor.  The 
motors  in  question  are  %  HP  machines  donated  by  the  Industrial  Electric  Shop  (lES)  on  the  E' 
District  Coast  Guard  base  in  Boston,  Massachusetts.  These  machines  were  previously  used  for 
intake  and  exhaust  ventilation  in  the  main  engine  room  aboard  one  of  the  Coast  Guard’s  49ft 
buoy  boats  (Eig.  3-2a).  Eigure  3-2b  shows  a  photograph  of  one  of  the  motors,  and  Table  3-1 
summarizes  key  motor  parameters. 


Figure  3-2:  (a)  One  of  the  Coast  Guard’s  49ft.  buoy  boats  and  (b)  one  of  its  exhaust  fans.  The  fans  nsed  in 
this  thesis  are  identical  to  the  one  pictured. 


Table  3-1:  Marine  Motor  Data 


Make:  Marathon  Electric  /  Marine  Duty  IEEE 

Model:  4VE  56C34  E5506A 

Type:  SCS 

Design:  B 

Code:  N 

Bearings:  6302 

Voltage:  120V 

Current:  10. 8  A 

RPM:  3450 

Single  phase 

60  Hz 

y4Hp 

#  Rotor  Slots:  34 

Poles:  2  (1  pair) 
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The  following  sections  describe  the  experimental  verification  of  each  of  the  key  parts  of 
the  process  described  above.  The  issues  relevant  to  this  thesis  are: 

•  Speed  estimation  using  rotor-slot  harmonics 

•  Detection  and  use  of  eccentricity-related  signals 

•  Detection  and  use  of  characteristic  bearing  frequencies 


3.5.1  Estimating  Speed  with  Slot  Harmonics 

An  accurate  speed  estimate  is  key  to  the  monitoring  process  described  in  this  chapter. 
Speed-related  current  modulations  arise  from  the  interaction  of  variations  in  the  air-gap 
permeance  with  the  air-gap  magnentomotive  force  (mmf).  Variations  in  the  air-gap  permeance 
are  caused  by  rotor  slotting  and  rotor  eccentricity.  It  can  be  shown  that  these  effects  result  in 
current  harmonics  at  the  following  frequencies  [24] : 

fsh  =  fi  ±  rid)  (^)  ±  ^w]  (3-1) 


In  Eq.  3-1,  /;  is  the  supply  frequency,  k  is  an  integer  in  the  range  {0,1,2...};  R  is  the 
number  of  rotor  slots;  Ud  =  0,+l,...,  is  the  order  of  rotor  eccentricity;  s  is  per  unit  slip;  p  is  the 
number  of  poles  pairs;  and  n^^,  =  +l,+3,...,  is  the  air-gap  mmf  harmonic  order  [24].  Given  values 
for/;,  7?,  and  p,  assuming  k  =  I  and  static  eccentricity  (nd=0),  one  can  use  these  slot  harmonics  to 
estimate  the  slip  and  subsequently  the  speed  from  the  equation: 


s  = 


Ns-N 

Ns 


(3-2) 


where  Ns  is  the  synchronous  speed^  of  the  machine  in  RPM  and  N  is  the  actual  speed  in  RPM. 
Initial  testing  of  the  speed  estimation  process  was  performed  using  the  test  machine  described 
above.  The  current  procedure  follows  three  steps: 


1.  Determine  a  range  of  values  for  four,  low-order  rotor-slot  harmonics  using  Eq.  3- 
1  and  n^^,  in  the  range  between  -3  and  3  (low  order).  An  initial  slip  estimate  is  also 
required.  In  our  initial  tests  this  estimation  is  obtained  from  direct  measurements, 
but  one  could  also  use  the  rated  speed  written  on  the  boiler  plate. 

2.  Examine  the  current  spectrum  in  the  range  calculated. 

3.  “Lock  on”  to  the  four  low  order  slot  harmonic  spectral  components. 

4.  Determine  speed  estimate. 


*  Synchronous  speed  is  the  number  of  pole  pairs  divided  by  3600,  (Ns  =  3600/p). 
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To  demonstrate  this  process,  several  experiments  were  performed  using  the  test  motor 
identified  in  Table  3-1.  In  each  experiment,  the  speed  of  the  motor  was  changed  by  placing  a 
different  weight  on  one  of  the  fan  blades.  Table  3-2  lists  the  values  needed  to  estimate  the 
locations  of  the  four  slot  harmonics  discussed  previously.  Given  an  initial  speed  estimate  of 
approximately  3500RPM,  the  four  slot  harmonics  are  expected  to  appear  in  the  range  between 
1800Hz  and  2200Hz.  This  is  clearly  indicated  in  Fig.  3-3. 


Table  3-2:  Slot  Harmonic  Variables 


Fi  =  60  Hz 

k  =  l 

R  =  34 

nd  =  0 

p  =  2 

tim  =  +1 ,  +3 

Frequency:  Hz 

Figure  3-3:  Current  Spectrum  for  the  test  motor  in  range  from  1800  Hz  to  2200Hz.  The  signals  indicated 

with  ovals  are  the  four  rotor-slot  harmonics. 

The  low  order  rotor  slot  harmonics  are  analyzed  more  closely  below.  The  harmonics  are 
presented  in  amplitude  order,  Hw  =  being  the  strongest  and  -3  the  weakest.  Manual  speed 
estimates  were  recorded  using  a  handheld  tachometer.  Each  labeled  plot  of  individual  slot 
harmonics  is  prefaced  by  a  corresponding  table.  Equation  3-1  is  used  to  calculate  where  that 
particular  harmonic  is  expected  in  the  current  spectrum.  The  actual  slot  harmonic  location  is  also 
recorded  for  comparison.  In  each  case,  the  slot  harmonic  estimate  differs  less  than  two  tenths  of 
a  percent  from  the  expected  value. 
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Table  3-3:  Slot  harmonic  Comparison  for  =  +1 


Speed 

Expected  fsh 

Actual  fsh  = 

Reference 

3512 

2050.1 

2047.3 

With  lOg  added 

3511 

2049.6 

2046 

With  13g  added 

3510 

2049 

2045.5 

45 
40 
35 
30 
-g  25 

I  20 

15 
10 
5 
0 

2044  2045  2046  2047  2048  2049 

Frequency:  Hz 

Figure  3-4:  Slot  Harmonic  (n„  =  +1) 


Table  3-4:  Slot  harmonic  Comparison  for  =  -1 


Speed 

Expected  fsh 

Actual  fsh  = 

Reference 

3512 

1930.1 

1927.3 

With  lOg  added 

3511 

1929.6 

1926 

With  13g  added 

3510 

1929 

1926.5 

Figure  3-5:  Slot  Harmonic  (n„  =  -1) 
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Table  3-5:  Slot  harmonic  Comparison  for  =  +3 


Speed 

Expected  fsh 

Actual  fsh  = 

Reference 

3512 

2170.1 

2167.3 

With  lOg  added 

3511 

2169.6 

2166 

With  13g  added 

3510 

2169 

2166.7 

2165  2165.5  2166  2166.5  2167  2167.5  2168  2168.5 

Frequency:  Hz 


Figure  3-6:  Slot  Harmonic  (n„  =  +3) 


Table  3-6:  Slot  harmonic  Comparison  for  n„  =  -3 


Speed 

Expected  fsh 

Actual  fsh  = 

Reference 

3512 

1810.1 

1807.3 

With  lOg  added 

3511 

1809.6 

1806 

With  13g  added 

3510 

1809 

1806.4 
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There  is  a  small  but  relatively  eonsistent  shift  between  the  caleulated  value  and  aetual 
value  for  each  set  of  slot  harmonics.  This  shift  could  be  due  to  a  slight  variance  in  supply 
frequency  ifi)  or  the  result  of  manual  speed  measurement.  This  particular  test  was  performed  at 
approximately  1:00pm  on  a  week  day,  a  time  when  the  power  grid  typically  experiences  a 
heavier  load.  It  is  possible  that  the  supply  frequency  may  have  dipped  slightly  during  this  test. 
For  example,  a  frequency  dip  of  0. 1  Hz  results  in  a  3  Hz  decrease  in  slot  harmonic  frequency  in 
the  current  spectrum.  In  order  to  develop  an  automated  speed  detection  algorithm  it  may  be 
necessary  to  determine  the  expected  supply  frequency  variation.  This  problem  could  be 
eliminated  altogether  by  using  the  current  spectrum  to  determine  supply  frequency  instead  of 
relying  on  user  input.  In  marine  applications  this  variation  can  be  significant  and  should  be 
taken  into  account. 

It  appears  that  each  of  the  variables  in  Eq.  3-1  can  be  assumed  constant  given  a  specific 
situation  with  the  exception  of  slip.  To  initially  find  the  slot  harmonics,  a  search  window  could 
be  set  to  start  at  the  point  of  zero  slip  and  moved  down  the  current  spectrum.  Once  the  spike  is 
detected  that  value  of  slip  used  to  position  the  search  window  would  be  recorded  and  used  in 
near  real  time  diagnostic  calculations. 


3.6  Detecting  Faults  using  Current  Spectral  Quantities 
3.6.1  Mechanical  Imbalance 

The  first  source  of  harmonics  injected  into  the  stator  current  of  an  induction  machine 
originates  from  cyclical  load  torque  variation  or  torque  modulation.  These  position-dependent 
changes  in  torque  are  common  in  electromechanical  loads.  The  frequencies  at  which  these 
components  are  produced  in  the  stator  current  are  predicted  by  Eq.  3-3.  In  this  equation,  Fi  is 
the  fundamental  supply  frequency,  s  is  the  slip,  p  is  the  number  of  poles,  and  m  is  an  integer 
value  [26]. 


load 


1  +m 


l-s 

p 

2 


(3-3) 


Observing  these  stator  current  harmonics  is  the  first  step  toward  confirming  the  NILM’s 
ability  to  reliably  detect  more  complex  load  and  bearing  problems.  To  verify  the  NIEM’s  ability 
to  detect  cyclical  load  torque  variation  the  axial  fan  used  to  produce  the  slot  harmonics  above 
will  be  used  to  verify  two  modes  of  imbalance  and  verify  their  presence  with  the  vibration 
spectrum. 

Vibration  monitoring  is  currently  the  primary  method  used  to  detect  the  presence  of  fault 
conditions  in  induction  machines.  The  mechanical  vibrations  are  associated  with  variations  in 
the  physical  air  gap  of  the  machine  and  stator  currents  are  generated  at  predictable  frequencies 
related  to  the  electrical  supply  and  vibration  frequencies  [18].  This  section  will  confirm  the 
value  of  current  monitoring  for  induction  machines  by  correlating  the  relationship  between 
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vibration  and  current  frequencies  caused  by  known  imposed  induction  machine  faults  (e.g. 
imbalance  and  bearing  faults). 

To  record  vibration,  a  medium  size  vibration  sensor  was  mounted  on  the  outside  of  the 
axial  fan  housing  as  depicted  in  Fig.  3-8.  This  orientation  was  selected  using  the  standard 
portable  vibration  recorder  (STS- 101)  operating  procedure  employed  by  the  Coast  Guard 
intermediate  maintenance  activity  [27].  This  sensor  was  not  used  to  provide  empirical  amplitude 
data  but  simply  used  to  establish  position  in  the  vibration  spectrum  and  show  relative  amplitude 
changes  between  known  conditions.  This  specific  “geo”  sensor  is  large  enough  to  produce  a 
current  that  can  be  wired  directly  into  a  NILM  channel  and  processed  in  parallel  with  current 
data.  Thus,  both  vibration  data  and  current  spectrum  data  will  be  presented  in  parallel  and 
represent  identical  moments  in  time. 


Figure  3-8:  Vibration  Sensor  Mount 


Equation  3-4  and  3-5  below  provide  a  framework  for  implementing  two  different  modes  of 
imbalance,  mass  variation  and  radial  location  variations. 

"^load  ~  "^Initial  T ^Sin(^Tlf  eccO  (^"4) 


Ti  =  mgr  (3-5) 

It  is  expected  that  an  increase  in  Tj  due  to  an  increase  in  mass  or  radial  distance  from  the 
axis  of  rotation  will  result  in  a  proportional  increase  in  Tioad-  Unbalanced  loads  will  cause  motor 
torque  to  vary  as  indicated  where  m  is  mass  of  the  imbalance,  g  is  the  gravitational  constant,  and 
r  is  the  radial  distance  of  the  mass  from  the  axis  of  rotation. 

This  first  test  looks  at  the  NILMs  ability  to  detect  cyclical  load  torque  variations  by 
varying  mass  and  keeping  the  location  of  the  mass  or  radius  constant.  Ti  is  expected  to  vary 
proportionally  with  mass  resulting  in  a  proportional  increase  in  the  total  torque  load  of  the  motor. 
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Figure  3-9  shows  the  methodology  used  to  imbalance  the  impellor.  Because  the  impellor  is 
spinning  at  nearly  3600  RPM  it  was  important  to  ensure  that  the  mass  was  secure.  Simply 
adding  bolts  of  varying  length  to  vacant  treaded  holes  in  the  hub  was  the  safest  and  most 
economical  way  to  perform  these  tests. 


Figure  3-9:  Axial  Fan  Imbalance  Method  (magnitude) 


Entering  Eq.  3-3  with  known  parameters,  supply  frequency  (/)  =  60  Hz),  harmonic 
number  {k  =  1),  number  of  poles  (p  =  2)  and  slip  (s)  is  calculated  with  Eq.  3-2  to  be  0.024, 
resulting  in  a  calculated /^cc  =  118.5  Hz.  The  relationship  between  physical  vibrations  and  the 
current  spectrum  can  be  expressed  by  Eq.  3-6  where  k  is  chosen  to  be  1. 


fecc  ~  \fl  i  k  *  fvibration\  (3-6) 

Entering  Eq.  3-6  with  known  values  for  fecc,  fi,  and  k,  the  frequency  of  vibration  can  be 
predicted  at  1 18.5  -  60  or  58.5  Hz  in  the  vibration  spectrum. 

Three  tests  were  run  in  this  set,  the  first  was  a  reference  (no  weight  added),  a  test  with  7 
grams  was  added  1”  away  from  the  axis  of  rotation,  and  a  test  with  10  grams  added  1”  away 
from  the  axis  of  rotation.  Eigure  3- 10a  and  b  presents  the  resulting  vibration  and  current 
spectrums  respectively.  This  test  confirms  the  expected  location  and  amplitude  relationship 
using  the  vibration  spectrum. 


50 


Vibration 


Blue  Current 


o'  "  *  '  - ^ - 1  ^  ^ 

116  117  118  119  120  121  122 

Frequency:  Hz 


Figure  3-lOa  and  b:  Varying  Mass  Test  Vibration  and  Current  Eccentricity 


The  next  test  performed  assess  the  NILM’s  ability  to  monitor  an  induetion  machine 
imbalance  using  a  constant  mass  but  varying  the  location  with  respect  to  the  center  of  axial 
rotation.  A  6.8  gram  calibrated  fan  weight  from  a  vibration  analysis  kit  at  the  Coast  Guard 
Industrial  Electrical  Shop  was  used  to  induce  a  vibration  as  shown  in  Fig.  3-11.  The  weights  are 
“C”  type  which  can  be  moved  along  the  leading  edge  of  the  impellor  blades.  Three  tests  were 
recorded,  the  first  test  represented  a  reference  with  no  weight  added,  the  second  placed  the 
weight  3”  from  the  center  of  rotation  and  the  final  test  moved  the  same  weight  to  a  position  3.25” 
from  the  center  of  rotation.  Figure  3-11  on  the  following  page  presents  a  weighted  run  with  the 
weight  3”  from  the  axis  of  rotation.  This  test  is  significant  because  this  showed  how  the  NIFM 
could  be  potentially  used  to  detect  misalignments  and  changes  in  alignment. 


Figure  3-11:  Fan  with  6.8  gram  weight  3  inches  from  center 
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Figure  3-12:  Axial  fan  test 


A  laser  tachometer  was  used  to  monitor  the  speed  of  the  fan  for  each  run.  The  speed  of 
the  fan  remained  relatively  constant  (3507  RPM)  over  each  test.  Using  Eq.  3-3,  the  eccentricity 
frequency  (fecc)  for  the  fan  was  calculated  to  be  118.45  Hz.  According  to  Eq.  3-2  and  3-3,  a 
linear  relationship  should  exist  between  the  eccentricity  component  amplitude  and  imbalance 
weight  distance  from  axis  of  rotation  or  radius  r. 


116  117  118  119  120  121  122 

Frequency:  Hz 

Figure  3-13:  Eccentricity  comparison  based  on  radial  location 

Eigure  3-13  presents  the  results  of  the  three  tests.  If  these  measured  eccentricities  are 
assumed  to  be  linear  then  the  total  system  response  can  be  considered  a  sum  of  simple  linear 
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responses.  This  test  eonfirms  the  NILM’s  ability  to  deteet  radial  ehanges  in  weight  imbalance 
and  it  reflects  the  proportionality  of  Eq.  3-4  and  3-5  within  reason. 


3.6.2  Motor  Bearing  Damage  Detection 

Numerous  studies  have  shown  that  bearing  problems  account  for  over  40%  of  all 
machine  failures  in  large  machinery  applications  [16],  [25].  There  is  strong  evidence  that  stator 
current  monitoring  can  provide  bearing  health  indications  similar  to  typical  vibration  analysis 
without  having  to  access  the  motor.  A  significant  amount  of  work  has  been  done  demonstrating 
the  feasibility  of  detecting  specific  bearing  faults  by  correlating  characteristic  bearing 
frequencies  to  the  spectral  components  of  the  stator  current  [18], [28].  Monitoring  of  roller 
bearings  is  very  important  in  applications  where  failure  would  result  in  loss  of  human  life  or 
replacement  would  require  long  production  downtimes. 

This  series  of  tests  seeks  to  confirm  that  it  is  possible  to  use  characteristic  bearing 
frequencies  in  the  current  spectrum  to  find  bearing  degradation.  This  is  possible  because 
defective  bearings  generate  mechanical  vibrations  at  the  rotational  speeds  of  each  component  in 
the  bearing.  These  characteristic  frequencies,  which  are  related  to  the  inner  or  outer  raceways 
and  /  or  balls  or  roller  geometry,  can  be  calculated.  Still  today  mechanical  vibration  analysis 
techniques  are  used  to  monitor  these  frequencies  in  order  to  determine  machine  bearing  health. 
The  relationship  of  the  bearing  vibration  to  the  stator  current  spectrum  can  be  determined  by 
recalling  that  any  air  gap  eccentricity  produces  anomalies  in  the  air  gap  flux  density  [29].  The 
test  in  this  section  attempt  to  confirm  that  the  NILM  can  detect  the  initial  signs  of  bearing 
damage  and  ultimately  bearing  failure. 

By  knowing  bearing  geometry  and  the  speed  of  the  induction  machine,  specific  bearing 
problems  can  be  predicted  and  located  (e.g.  inner  /  outer  raceway,  specific  bearings,  etc.).  The 
same  vent  fan  used  in  the  last  test  was  monitored  for  bearing  wear  testing. 


3. 6.2.1  Inducing  Realistic  Bearing.  Faults 

There  are  a  myriad  of  conditions  that  cause  bearing  faults  in  industry.  Some  of  the  most 
common  faults  are  flaked  running  surfaces  due  to  material  fatigue,  excessive  clearance  due  to 
extensive  wear,  excessive  preload  conditions,  and  poor  lubrication.  Bearing  material  fatigue 
usually  results  in  a  crack  created  in  rings  and  rolling  elements  and  originates  from  a  pre-damaged 
spot  on  the  raceway  /  bearing  or  because  of  an  inhomogeneous  cast  during  manufacture. 
Damage  starts  below  the  surface  and  propagates  to  the  surface  of  the  rolling  element  where  it 
causes  material  to  break  loose  forming  pitting  in  the  contact  areas  as  shown  in  Fig.  3-14  [29]. 
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3-14:  Example  of  Outer  Race  Spalling  in  a  Deep  Groove  Ball  Bearing 


3.6.2.2  Realistically  Damaging  the  Inner  and  Outer  Bearing  Race 

Most  small  to  medium  bearing  assemblies  are  machine  assembled  and  sealed  to  preserve 
the  manufacturer’s  lubricant  and  prevent  foreign  object  intrusion.  This  makes  it  a  challenge  to 
damage  a  specific  location  in  a  bearing  race  without  affecting  other  bearing  elements.  A  custom 
clamp,  vice  and  Dremel  with  a  diamond  tip  bit  was  used  to  damage  the  inner  and  outer  race  in 
two  separate  6203  bearings  in  a  manner  very  similar  to  that  seen  in  Fig.  3-14.  Figure  3- 15a,  b, 
and  c  are  photos  of  the  damaging  process. 


Figure  3-15a,b,c:  Bearing  Damaging  Process 

The  bearings  were  held  in  place  using  a  hand  vice  and  stationary  vice  while  a  small 
diamond  tip  bit  was  negotiated  around  the  raised  race  edges  and  bearing  cage.  The  bit  was  held 
in  place  long  enough  to  remove  the  hardened  surface  of  the  race  to  a  depth  of  between  1  and  2 
millimeters.  After  damaging  the  inner  and  outer  races  the  bearing  need  to  be  operated  to  allow 
time  for  fault  development.  Each  of  the  damaged  bearings  was  initially  new  with  tight  bearing 
cages  and  race  fits.  In  order  to  wear  two  sets  of  damaged  bearings  simultaneously,  the  damaged 
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inner  race  bearing  was  installed  in  the  non-load  end  and  the  damaged  outer  race  bearing  was 
installed  in  the  shaft  end  of  the  axial  fan  motor.  The  fan  was  energized  for  approximately  10 
days.  Next,  the  motor  was  disassembled  and  the  damaged  bearings  installed  on  the  non-load  end 
(bearings  with  inner  race  damage)  were  replaced  with  a  new  set  and  the  motor  reassembled. 


Figure  3-16:  Axial  Motor  Test  Stand 


The  NILM  was  configured  to  monitor  six  channels,  two  recording  vibration  data  like  that 
used  in  section  3.2.2  as  well  as  voltage  and  current  data  for  two  motors. 

3.6.23  Monitoring  Bearing  Race  Damage 

A  derivation  for  the  bearing  vibration  equations  can  be  found  in  [29].  The  solution  to  the 
inner  and  outer  race  bearing  frequency  components  can  be  found  in  the  vibration  spectrum 
where  n  is  the  number  of  balls  in  the  bearing,  PD  is  the  bearing  pitch  diameter,  BD  is  the  ball 
diameter,  (3  is  the  contact  angle  of  the  balls  on  the  races  and  /„„  is  the  mechanical  rotor  speed  in 
Hertz.  Pitch  diameter,  ball  diameter,  contact  angle  and  the  number  of  bearings  can  be  attained 
by  the  manufacturer.  Once  determined,  this  value  must  be  entered  into  Eq.  3-9  to  determine 
current  spectrum  location. 


(3-7) 


(3-8) 
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The  manufacturer  specifics  for  the  6302  bearing  assembly  installed  are  listed  below: 

•  PD  =  28.5  mm 

•  BD  =  7.5  mm 

•  n  =  8 

•  =  0  degrees* 

*  After  inducing  damage,  bearings  become  looser.  Given  the  direction  of  force  due  to 
axial  fan,  it  appears  likely  that  p  is  between  10  and  15  degrees  in  these  bearings. 


The  fan  speed  varied  from  3510  -  3520  RPM  with  the  damaged  inner  and  outer  race 
bearing  test.  It  is  important  to  note  that  the  fans  were  oriented  horizontally  for  the  following 
outer  and  inner  bearing  race  damage  test.  Given  machine  speed,  the  frequency  of  bearing 
vibration  was  calculated  at  173Hz  and  the  corresponding  current  spectrum  harmonics  for  outer 
race  faults  was  calculated  at  113  Hz  and  233  Hz  for  k  =  and  k  =  +1  respectively  using  Eq.  3-9. 

fbrg^\fl±k*  fy\  (3-9) 


Vibration 


Figure  3-17:  Outer  Race  Frequency  Vibration  Spectrum 

The  Figure  3-17  shows  the  vibration  spectrum  in  the  vicinity  of  the  calculated  outer  race 
frequency.  Figure  3-18,  shows  the  current  spectrum  in  the  vicinity  of  the  sum  and  difference  of 
60  Hz  and  fy.  It  is  difficult  to  determine  if  a  corresponding  signature  could  be  identified  in  the 
current  spectrum  in  either  plot. 
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Motor  Current 


Motor  Current 


Motor  Current 


Figure  3-19:  Current  Spectrum  (230  -  242Hz) 


Unfortunately,  the  eharacteristic  outer  race  resides  very  close  to  the  second  harmonic  of 
the  eccentricity  frequency  used  earlier  in  this  chapter  to  detect  machine  imbalance.  This  is  a 
problem  that  could  be  remedied  if  considered  in  the  system  design  phase. 

The  same  process  was  used  to  present  the  damaged  inner  race  data.  Figure  3-20  and  21 
on  the  following  page  presents  a  comparison  of  a  good  bearing  and  a  bearing  with  a  damaged 
inner  ring  located  on  the  load  end  of  the  induction  machine.  Referring  back  to  Eq.  3-7,  the  inner 
race  characteristic  frequencies  are  calculated  to  be  235  Hz  and  355  Hz  for  k  =  -I  and  k  =  +\ 
respectively.  The  inner  race  damage  showed  up  much  better  than  the  outer  race  data  and  is  in  an 
optimal  location  in  the  current  spectrum. 
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Motor  Current 


Figure  3-20:  Inner  Bearing  Race  Spectral  Component  (k  =  -1) 


Motor  Current 


Figure  3-21:  Inner  Bearing  Race  Spectral  Component  (k  =  +1) 

It  is  very  difficult  to  damaging  bearings  in  a  precise  manner.  The  difference  between  the 
amplitudes  of  the  inner  and  outer  race  spectral  components  may  be  attributed  to  a  number  of 
factors  such  as  system  or  equipment  resonances.  Regardless,  the  inner  ring  data  appears  strong 
compared  to  spectral  comonents  in  the  vicinity  of  each  inner  ring  bearing  frequency.  In  each 
presented  case,  the  damaged  bearing  was  installed  on  the  shaft  end  of  the  rotor. 
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3.7  Impediments  to  Detecting  Faults 

3.7.1  Pressure 

A  considerable  number  of  lessons  were  learned  while  performing  the  test  in  this  ehapter. 
The  first  involved  test  in  section  3.6.1,  mechanical  imbalance.  Originally  the  goal  was  to  load 
the  motors  by  restricting  the  airflow  into  the  motor  under  test.  Results  were  found  to  be 
ineonelusive  because  this  restriction  caused  such  a  significant  change  in  pressure  across  the  fan 
(difference  between  pressure  at  the  inlet  and  atmospheric  pressure  in  this  case)  that  the 
eeeentrieities  were  severely  damped.  In  comparing  Fig.  3-22  on  the  following  page  and  Fig.  3- 
10b,  amplitude  and  separation  is  significantly  affected  by  amount  of  pressure  across  an  axial  fan. 


Current  Spectrum  Blue 


Frequency:  Hz 

Figure  3-22:  Pressure  Damped  Eccentricity 


3.7.2  Orientation 

Gravity  plays  a  significant  role  in  deteeting  spectral  quantities  that  are  a  funetion  of 
torque  variation  as  well  as  eccentricity.  This  is  evident  considering  Eq.  3-5,  but  the  not  intuitive 
as  to  how  dramatie  the  effect  might  be.  Initially,  the  bearing  damage  tests  were  conducted  with 
the  fans  oriented  vertically.  Each  test  was  repeated  several  times  and  a  week  of  operation  added 
to  the  bearing  sets  to  develop  the  damage.  Time  and  deepening  of  the  race  damage  did  little  to 
affect  the  results.  Eigure  3-23  on  the  following  page  shows  very  different  results  while  using  the 
same  bearings  installed  for  the  Eig  3-20  and  3-21  tests.  The  eharaeteristic  bearing  frequencies  at 
231.5  Hz  and  351.5  are  hardly  deteetible.  In  this  test  the  bearings  were  installed  in  the  shaft  end 
of  the  rotor  and  the  rotor  positioned  on  top. 
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Figure  3-23:  Damaged  Inner  Race  Bearing 
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3.8  Conclusion 

Multiple  modes  of  motor  damage  ean  be  deteeted  using  non-intrusive  techniques.  There 
are  a  significant  number  of  factors  that  will  contribute  to  the  successful  implementation  of  non- 
intrusive  monitoring  techniques.  If  current  spectral  quantities  are  to  be  used  for  diagnostics, 
effort  must  be  made  to  manage  the  relationship  between  physical  design  characteristics  and 
harmonics  of  60Hz.  Using  the  example  of  the  axial  fan  above,  the  inner  bearing  race 
characteristic  frequency  (235Hz)  was  very  close  to  the  fourth  harmonic  of  power  (240Hz).  This 
would  make  automating  the  process  depicted  above  more  difficult.  Every  effort  should  be  made 
to  take  these  types  of  design  decisions  into  account. 

Machine  loading  significantly  affects  the  signature  produced  by  eccentric  spectral 
quantities  in  the  current  spectrum.  Again,  automating  a  process  to  monitor  and  trend  these 
quantities  becomes  much  more  difficult  if  the  amplitude  of  these  quantities  varies  significantly. 

Induction  machines  are  just  one  example  of  a  system  that  current  spectral  quantities  could 
be  used  to  develop  machine  health  diagnostics.  Monitoring  techniques  are  continually  being 
improved.  It  is  for  this  reason  that  current  monitoring  should  remain  a  source  of  diagnostic 
information. 


60 


4  Multiple  load  monitoring 

4.1  Laboratory  Testing 


The  initial  goal  in  developing  the  proeess  of  using  a  single  aggregate  sensor  to  monitor 
multiple  loads  was  to  start  simple.  Two  axial  fans  like  those  used  in  ehapter  3  with  the  same 
make  and  model  were  arranged  as  shown  in  Fig.  4-1.  The  test  stand  was  built  with  using 
standard  industrial  shelving,  additional  %  ineh  plywood  for  additional  support,  and  HVAC  foam 
board  to  enclose  part  of  the  structure.  A  local  MIT  HVAC  post  doctoral  researcher,  Peter 
Armstrong,  was  consulted  regarding  ventilation  intake  and  exhaust  motor  configuration  and 
loading. 


Figure  4-1:  Multi-motor  Test  Stand 

United  States  Coast  Guard  Technical  Publication  2830B  (WMEC  905  Class  Ventilation 
Systems  Manual)  was  then  referenced  to  apply  a  main  engine  room  intake  and  exhaust 
relationship  to  a  small  scale  lab  test.  In  most  cases  the  main  space  intake  and  exhaust  fans  were 
located  in  the  space  they  serviced.  Therefore,  the  intake  fan  experienced  the  drag  of  ducting  on 
the  intake  side  of  the  individual  motor.  Conversely,  the  exhaust  fans  typically  experienced 
slightly  less  drag  on  the  discharge  side  of  the  exhaust  fan.  This  is  significant  because  the 
difference  in  loading  provides  an  opportunity  in  an  industrial  setting.  Many  times  it  is  cheapest 
to  use  the  same  motor  for  both  intake  and  exhaust  in  industrial  or  marine  applications  because  of 
maintenance  and  supply  chain  concerns.  In  this  case,  two  identical  motors  loaded  differently 
will  result  in  each  motor  having  a  characteristic  slip  that  can  be  used  to  identify  each  from  a 
single  aggregate  source. 
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It  was  determined  that  the  intake  fan  should  be  configured  to  experience  a  flow  resistance 
similar  to  that  resulting  from  a  half  blocked  motor  inlet.  The  exhaust  fan  should  experience  the 
flow  resistance  resulting  from  a  quarter  blocked  discharge  to  simulate  the  relationship  described 
above.  Figure  4-2  is  a  simple  diagram  of  the  resulting  setup.  The  material  located  above  the 
vent  fans  was  installed  to  prevent  flow  oscillations. 


Figure  4-2:  Initial  Multi-motor  test-stand 


4.1.1  NILM  Setup 

The  NILM  setup  in  Fig.  4-3  was  configured  to  record  six  channels  of  data.  The  first  two 
channels  were  used  to  record  vibrations  of  the  exhaust  motor  for  the  reasons  expressed  in  chapter 
2.  Channels  3  and  4  were  used  to  record  aggregate  voltage  and  current  respectively  and  channels 
5  and  6  recorded  only  the  blue  motor.  A  block  diagram  of  this  test  setup  is  shown  in  Fig.  4-4. 


Figure  4-3:  NILM  Setup  for  Aggregate  Sensor  Testing 
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Figure  4-4:  Block  Diagram  for  Multi-motor  NILM  Test 

The  black  node  represents  a  voltage  connection,  an  unfilled  circle  represents  a  current 
transducer  and  the  highlighted  circle  on  the  blue  motor  represents  the  vibration  sensor.  It  was 
necessary  to  use  channels  3  and  5  for  voltage  in  order  to  prevent  a  significant  time  lag  between 
either  aggregate  or  blue  motor  voltage  current  pairs.  This  is  important  for  preprocessing  data  in 
the  future. 


4. 1 .2  Modifying  the  Single  Motor  Process 

The  primary  difference  between  the  single  load  block  diagram  in  Fig.  3-1  and  Fig.  4-5  is 
that  key  data  pertaining  to  each  individual  load  (e.g.  number  of  rotor  slots,  number  of  pole  pairs, 
etc.)  is  organized  by  load  into  arrays.  These  parameters  are  used  to  identify  slot  harmonics 
associated  with  a  specific  machine,  calculate  slip  for  that  machine  and  pass  that  information  with 
all  other  information  known  about  that  particular  machine  to  the  diagnostics  section.  Recall  in 
the  previous  example,  the  GUI  would  have  a  data  entry  function  where  the  key  motor  parameters 
for  each  machine  under  test.  This  parametric  data  is  fed  into  the  Information  and  Computation 
function  and  organized  according  to  the  specific  induction  machines  feed  by  the  load  center 
being  monitored.  This  function  stores  the  key  parameters  and  uses  the  slot  harmonic  calculation 
to  determine  slip  for  each  induction  machine.  It  then  packages  this  information  in  vectors  and 
passes  it  to  the  diagnostic  subroutines  according  to  preferences  entered  into  the  GUI.  This 
process  is  repeated  until  each  of  the  entered  loads  is  initially  established  and  passes  through  the 
specified  diagnostics.  Figure  4-5  on  the  following  page  presents  the  multi-motor  approach  using 
a  general  block  diagram.  There  are  a  number  of  advantages  to  using  such  a  process.  First,  the 
graphical  user  interface  (GUI)  could  be  used  to  enter  the  number  of  loads,  type  of  loads  and  key 
attributes  of  those  loads  into  a  computational  module.  Additionally,  the  GUI  could  also  be  used 
to  establish  or  eliminate  diagnostic  test  to  performed  on  a  specific  machine  or  load. 
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Figure  4-5:  General  Multi-motor  Diagnostic  Process  Block  Diagram 


4. 1 .3  How  test  were  performed 

The  PCI  1710  was  set  to  sample  at  a  rate  of  48.456  kHz  or  8.076  kHz  per  channel  for  90 
seconds.  Once  the  PCI  1710  card  was  configured  and  destination  directory  set,  both  motors  were 
energized  and  data  recording  started.  After  45  seconds  the  red  motor  was  secured  leaving  only 
the  blue  motor  running.  At  each  of  the  two  steps,  current,  real  power,  reactive  power,  power 
factor,  axial  vibration  and  each  fan’s  RPM  were  recorded.  This  manual  data  was  recorded  and 
used  to  confirm  NILM  operation. 


4.2  Determining  Speed  of  Multiple  Machines  using  Aggregate  Sensor  data 

Ninety  seconds  of  raw  current  data  was  data  was  recorded  and  stored  for  each  test 
conducted.  To  determine  the  feasibility  of  monitoring  these  two  induction  motors  from  an 
aggregate  point,  the  eccentricity  test  procedure  from  Chapter  3  are  used.  Figure  4-7  shows  the 
recorded  current  in  the  time  domain.  The  window  in  the  figure  highlights  the  section  of  data 


.Aggregate  Ourrent 


Figure  4-6:  Multi-Motor  Time  Domain  Current 
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selected  to  create  the  frequency  domain  plots  or  spectral  current.  There  are  three  tests 
represented,  a  reference  case  with  no  imbalance,  a  case  with  10  grams  of  weight  added  one  inch 
from  the  axis,  and  a  run  with  13  grams  added  one  inch  from  the  axis.  The  time  domain  plot 
appears  all  red  because  the  larger  imbalance  draws  the  most  current.  Each  motor  is  virtually 
identical  with  the  exception  of  slip,  which  should  result  in  two  definitive  sets  of  slot  harmonics 
near  the  predicted  values.  Each  set  will  be  used  to  predict  the  speed  of  each  machine. 


The  same  methodology  discussed  in  Chapter  3  is  employed.  Recall  Eq  3-1  shown  below, 
the  slot  harmonic  equation  With  values  for  R,  and  p  known,  a  k  value  of  1  was  assumed  as 
well  as  static  eccentricity  (nd=0).  The  same  initial  values  in  Table  3-2  were  used. 


fsh  =  fi 


(kR  ±  rid) 


(^) 


± 


Table  4-1:  Slot  Harmonic  Variables 


Fi  =  60  Hz 

k  =  l 

R  =  34 

=  0 

II 

a, 

=■  +1?  ib3 

At  this  point  each  of  the  four  slot  harmonics  will  be  compared.  Tables  4-2  thru  4-5  list 
the  expected  location  and  actual  location  for  nw  =  -i-l,  -1,  -i-3,  and  -3,  respectively.  Both  motors 
corresponding  slot  harmonics  can  be  seen  at  predicted  locations  in  the  current  spectrum. 


Aggregate 
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4.2.1  Multi-load  Slot  Harmonic  Analysis 


The  tests  procedure  used  in  Chapter  3  was  performed  and  the  aggregate  current 
transducer  data  was  used  to  create  the  figures  below.  In  the  Uw  case,  the  current  transducer 
isolating  the  blue  motor  was  used  to  confirm  similar  amplitude  scheme  as  well  as  location  in  the 
current  spectrum.  They  do  in  fact  agree  as  do  the  remaining  three  harmonics. 


Table  4-2:  Slot  harmonic  Comparison  for  n„  =  +1 


Speed 

(Red) 

Speed 

(Blue) 

Expected 
fsh  (Red) 

Actual  fsh 
=  (Red) 

Expected 
fsh  (Blue) 

Actual  fsh  = 
(Red) 

Reference 

3482 

3508 

2033 

2032.5 

2047 

2047.4 

With  lOg  added 

3480 

3503 

2032 

2032.6 

2045 

2046 

With  13g  added 

3485 

3505 

2035 

2034 

2046 

2046.8 

Blue  Current  Transducer 


Frequency:  Hz 

Figure  4-8:  Slot  Harmonic  (n^  =  +1) 
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Table  4-3:  Slot  harmonic  Comparison  for  n„  =  -1 


Speed 

(Red) 

Speed 

(Blue) 

Expected 
fsh  (Red) 

Expected 
fsh  (Blue) 

Actual  fsh 
=  (Red) 

Actual  fsh  = 
(Blue) 

Reference 

3482 

3508 

1913 

1928 

1913 

1927.3 

With  lOg  added 

3480 

3503 

1912 

1925 

1913 

1926 

With  13g  added 

3485 

3505 

1915 

1926 

1914 

1927 

Aggregate 


Table  4-4:  Slot  harmonic  Comparison  for  n„  =  +3 


Speed 

(Red) 

Speed 

(Blue) 

Expected 
fsh  (Red) 

Expected 
fsh  (Blue) 

Actual  fsh 
=  (Red) 

Actual  fsh  = 
(Red) 

Reference 

3482 

3508 

2153 

2168 

2153.2 

2167.3 

With  lOg  added 

3480 

3503 

2152 

2165 

2153 

2166 

With  13g  added 

3485 

3505 

2155 

2166 

2154.2 

2166.8 

Aggregate 
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Table  4-5:  Slot  harmonic  Comparison  for  n„  =  -3 


Speed 

(Red) 

Speed 

(Blue) 

Expected 
fsh  (Red) 

Expected 
fsh  (Blue) 

Actual  fsh 
=  (Red) 

Actual  fsh  = 
(Blue) 

Reference 

3482 

3508 

1793 

1808 

1793.2 

1807.2 

With  lOg  added 

3480 

3503 

1792 

1805 

1793 

1806 

With  13g  added 

3485 

3505 

1795 

1806 

1794 

1806.5 

Aggregate 


Frequency:  Hz 

Figure  4-11:  Slot  Harmonic  (n„  =  -3) 

The  slot  harmonics  for  these  four  cases  show  up  relatively  consistently  within  a  Hertz  of 
the  expected  value.  Tachometer  readings  of  both  motors  varied  by  +  1  RPM.  Therefore,  an 
average  value  was  used.  The  most  significant  error  was  less  than  1/10*  of  a  percent,  which  is 
extremely  accurate  by  most  standards.  The  amplitude  does  not  seem  to  provide  any  definitive 
information  for  this  particular  set  of  test.  Additionally,  there  does  not  appear  to  be  a  relationship 
between  imbalance  and  slip. 

There  are  limitations  to  this  particular  monitoring  scheme.  It  is  proposed  in  this  chapter 
that  slot  harmonics  be  used  to  determine  speeds  of  multiple  motors  to  populate  a  slip  array  (si, 
S2,...,  s„)  where  n  =  number  of  motors  supplied  by  the  load  center.  In  this  scenario,  two  near 
identical  motors  are  operating  at  close  to  the  same  speed.  Consider  the  case  above  with  a  13 
gram  imbalance;  the  red  motor  is  defined  as  motor  1  and  blue  motor  defined  as  motor  2.  Using 
the  slot  harmonic  determined  slip  value  for  each  would  result  in  Si  =  0.9675  and  S2  =  0.9739. 
Inserting  those  values  into  the  cyclical  load  torque  equation  (Eq.  3-3)  used  earlier  and 
represented  below,  result  in  Fioad.i  =  1 18.05  and  fioad_2  =  1 18.4. 


load 


=  Fi 


1  +m 


V 

2 
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The  separation  between  these  two  spectral  components  is  less  than  one  half  of  a  Hertz 
and  would  be  nearly  impossible  to  reliably  differentiate  with  diagnostic  software.  Figure  4-13 
presents  the  actual  plot  of  the  aggregate  cyclical  load  component  scenario  described  above. 


Aggregate 


Frequency:  Hz 


Figure  4-12:  Scenario  Aggregate  Current  Spectrum  plot  of  Cyclical  Load  Component 

In  order  to  monitor  two  identical  motors  from  the  same  source,  it  will  be  necessary  to 
incorporate  a  subsequent  test  procedure  to  identify  which  of  the  two  identical  loads  was 
experiencing  irregularities.  The  NILM  could  still  alert  an  operator  and  /  or  technician  of 
irregularities  but  in  this  case  would  need  to  take  additional  actions  to  pinpoint  the  irregularity. 
The  same  can  be  shown  for  the  bearing  scenario  as  well. 

4.2.2  Design  Options 

This  test  using  two  near  identical  induction  machines  represents  a  very  difficult  scenario. 
However,  with  the  knowledge  of  how  the  current  spectrum  and  slot  harmonics  can  be  used  to 
identify  motors,  a  designer  could  choose  a  design  that  took  NILM  capabilities  into  account. 
Some  shipboard  HVAC  systems  are  already  moving  toward  more  NILM  friendly  architectures. 
Because  the  intake  and  exhaust  vent  fans  are  typically  loaded  differently,  some  naval  architects 
are  using  that  fact  to  choose  optimal  sizes  for  these  vent  fans.  The  majority  of  these  cases  results 
in  larger  intake  and  smaller  exhaust  vent  fans,  which  potentially  presents  a  better  opportunity  to 
differentiate  the  machines  using  the  current  spectrum. 

Many  new  ships  are  automating  main  engineering  space  HVAC  systems  altogether.  For 
example,  the  87  ft  Marine  Protector  Class  pictured  in  Fig.  4-14,  employs  a  24  inch,  2  Vi  HP  axial 
fan  intake  to  service  the  main  engine  room  and  a  15  inch,  1  HP  axial  exhaust  fan,  both  which  can 
be  seen  in  figure  4-15.  Each  is  operated  by  a  variable  speed  drive  which  can  be  set  manually 
with  the  panel  in  Fig.  4-16  or  automatically  based  on  main  engine  loading.  In  this  case,  it  is 
much  more  likely  that  these  fans  do  not  have  an  identical  number  of  rotor  slots  and  pole  pairs. 
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Figure  4-13:  USCGC  Flying  Fish 


Figure  4-14:  Fan  Control  Interface  Panel 


Figure  4-15:  Main  Engine  Room  Intake  and  Exhaust  Ventilation  Fans 
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4.3  Conclusions 


It  is  possible  to  use  speetral  quantities  to  monitor  multiple  induetion  maehine  loads  with 
one  eurrent  sensor.  Speeial  teehniques  may  be  required  in  order  to  determine  the  speeifie 
induetion  maehine  eausing  elevated  speetral  quantities  assoeiated  with  system  health  problems. 
The  key  faetors  in  monitoring  multiple  induetion  maehines  with  an  aggregate  sensor  are  number 
of  rotor  slots,  number  of  pole  pairs,  and  load  sehedule  or  environment  of  each.  In  an  ideal  world, 
a  diagnostic  for  monitoring  multiple  motors  from  one  load  center  would  be  able  to  monitor  a 
specific  cyclical  load  torque  spectral  quantities  for  each  induction  machine.  Depending  on  what 
the  induction  machine  is  driving,  this  may  or  may  not  be  possible. 
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5  Moving  to  a  new  NILM  Architecture 

5.1  NILM  System  Architecture 


This  chapter  will  analyze  the  general  NILM  arehiteeture  prior  to  June  2007,  when  this 
thesis  was  started,  and  present  a  new  arehiteeture  using  tools  and  knowledge  presented  in  the 
System  Design  and  Management  (SDM)  eurriculum.  The  original  arehiteeture  will  be  ealled 
first  generation  henceforth  and  the  emerging  arehiteeture  will  be  referred  to  as  seeond 
generation.  The  seeond  generation  architecture  presented  in  this  ehapter  represents  the  work  of 
MIT  Prof.  Steve  Leeb,  UNC  Prof.  Robert  Cox,  Dr.  Jon  Rodriguez  of  NemoMetries,  LCDR  Ethan 
Proper,  LT  Perry  Braneh,  and  the  author.  Refer  to  [19]  for  additional  teehnieal  details  relating  to 
the  first  generation  NILM.  The  resulting  arehiteeture  will  be  presented  at  the  same  teehnieal 
level  as  the  first  generation  NILM  is  introdueed. 

In  order  to  move  to  a  new  arehiteeture,  the  first  generation  hardware  will  be  analyzed  and 
deeomposed.  The  goal  of  the  NILM  system  from  a  eommereial  standpoint  is  to  be  developed 
into  a  robust  produet  that  delivers  value  at  a  reasonable  cost  to  a  wide  range  of  eustomers.  The 
NILM  system  has  the  potential  to  deliver  valuable  information  to  customers  in  numerous 
industrial  settings.  Beeause  of  the  variance  in  these  settings,  there  will  be  times  when  the 
priorities  must  be  balaneed,  eustomer  needs  and  manufacturability  for  example. 

5.1.1  Modeling  the  First  Generation  Architecture 

The  first  step  taken  in  analyzing  the  first  generation  arehiteeture  was  to  map  the  existing 
eoneept  using  a  general  modularized  0PM  or  Objeet  Proeess  Model.  The  intent  of  this  initial 
0PM  is  not  to  synthesize  the  arehiteeture  from  eoneept  but  to  show  first  and  seeond  level 
deeomposition  relationships. 

The  0PM  in  figure  5-1  was  developed  using  details  from  [19]  and  [20].  The  first 
generation  NILM  setup  uses  the  voltage  measuring  board  and  eurrent  transdueer  together  with  a 
power  supply  and  wiring  terminal  board  inside  a  NEMA  enelosure.  Outputs  from  the  voltage 
board  and  eurrent  transdueer  are  eonneeted  to  the  wiring  terminal  board.  The  terminal  board 
uses  a  Small  Computer  System  Interfaee  (SCSI)  or  “skuzzy”  port  to  interfaee  with  the  PC.  This 
port  is  eonneeted  direetly  to  the  analog  to  digital  eonverter  in  the  PC.  A  PCI  1710  analog-to- 
digital  data  eard  used  eonverts  an  analog  signal  in  the  range  of  -5  V  to  -i-5V  to  digital  "eounts"  in 
the  range  of  0  to  4096  amplitude  levels.  The  inside  of  a  standard  first  generation  NILM  NEMA 
enelosure  ean  be  seen  in  Eig.  5-1.  The  PC  is  a  standard  desktop  maehine  using  a  standard 
Pentium  proeessor,  RAM,  CD  ROM,  LCD  display  and  keyboard.  A  eomplete  system  installed 
aboard  the  CGC  Escanaba  ean  be  seen  in  seetion  5.1.2. 
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Transducers 


Terminal  Block 


Figure  5-1:  First  Generation  NILM  NEMA  Enclosure 


Current 
transducers 
conductors  to 
NEMA 


SCSI  port  to  PC 


Current  transducers  are  located  in  the  host  or  monitored  load  center  and  fed  into  the 
NEMA  enclosure  and  ultimately  to  the  terminal  block. 


Eigure  5-2  represents  the  first  order  0PM  that  depicts  these  first  and  second  level 
relationships. 


Figure  5-2:  Modular  Relation  OPM 
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Stepping  through  the  0PM  development  for  the  first  generation  NILM  system  resulted  in 
a  better  understanding  of  sub-system  boundaries.  For  example,  it  became  clear  that  a  boundary 
would  need  to  be  drawn  between  analog  and  digital  data  sets  even  though  they  could  have  been 
generally  referred  to  as  “data.” 


5.1.2  NILM  System  Decomposition 

In  order  to  enter  a  design  tool  such  as  a  Dependency  Structure  Matrix  a  careful 
decomposition  of  the  first  generation  system  must  be  performed.  Figure  5-3  shows  a  standard  NILM 
installation.  The  current  complete  NILM  system  is  comprised  of  commercially  available 
hardware  that  is  relatively  easy  to  acquire,  assemble,  and  install.  Commercial  or  COTs  hardware 
was  initially  chosen  because  of  its  wide  availability  and  low  cost.  A  standard  first  generation 
NILM  system  consists  of  a  NEMA-type  enclosure  to  house  the  voltage  transducers  and  their 
associated  power  supply,  a  personal  computer  (PC)  with  keyboard  and  monitor,  and  an 
uninterruptible  power  supply.  The  personal  computer  used  in  NILM  setups  runs  a  Linux 
operating  system.  Linux  is  a  very  reliable  system,  which  is  especially  important  for  system 
robustness.  Installation  of  a  NILM  setup  can  be  accomplished  using  readily  available  hardware 
and  tools.  For  installations  discussed  in  this  thesis,  the  hardware  was  attached  to  panels, 
bulkheads,  or  mounting  plates  with  nuts  and  bolts. 


NILM  Box  / 
Kev  board 


9 


9 


9 


9 


PC/CPU 


Display 


CHT 


Load 


Panel 


5-3:  Typical  NILM  Equipment  Installation  (CGC  Escanaba) 
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The  personal  eomputer  (PC),  display,  uninterruptable  power  supply  (UPS),  and  key  board 
are  all  standard  COTS  hardware.  The  NILM  box  however  has  not  had  the  same  level  of  produet 
development.  Moving  to  a  new  architecture  provides  the  team  the  opportunity  to  take  a  fresh 
look  at  customer  needs.  A  block  diagram  of  the  NILM  unique  hardware  is  presented  in  Fig.  5-4. 


Ccrrvnt  Cutpjt 


Fpow  S^io  S«p*vice  Power 

5-4:  Current  NILM  Block  Diagram 

A  modern  Design  Structure  Matrix  or  DSM  will  be  used  to  suggest  a  new  architecture 
and  will  require  a  preliminary  decomposition.  The  decomposition  will  need  to  dip  down  two 
levels  prior  to  entering  the  DSM.  Given  that  level  0  is  the  NILM  system  itself,  the  level  1  and  2 
decomposition  can  be  defined  as  depicted  in  Table  5-1. 


Table  5-1:  First  Generation  NILM  Decomposition 


First  Generation  NILM  Levei 

0  Decomposition 

Level  1  Decomposition 

Level  2  Decomposition 

System 

Pentium  Based  PC 

Power  SuddIv  (120  -  12V) 

CPU 

RAM 

CD  ROM 

Hard  Disk 

Video  Card 

Analog  -  Digital  Card 

Software 

Enclosure 

Display 

Power  SuddIv  (120  -  12V) 

Liauid  Crystal  DisDiav 

Enclosure 

UPS 

Battery 

Battery  Charging  System 

CPU  /  Converter 

Enclosure 

NILM  (NEMA  Enclosed  hardware) 

Power  SuddIv  (120  -  +  12  V) 

Voltage  Board 

Terminal  Board 

Enclosure 

Current  Transducers 

Current  Transducers 

Standard  Outlet  * 

Single  Phase  120  V  Power 

Monitored  System  * 

Voltage  Interface 

Current  Interface 

Enclosure 

*  Only  Applicable  parts  of  these  elements  listed 
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The  goal  in  defining  specifie  level  decomposition  is  to  keep  the  parts  at  approximately 
the  same  level  of  complication.  It  is  important  to  note  that  interconnections  between  the  parts 
are  not  listed  in  the  decomposition.  Standard  and  nonstandard  subsystem  interfaces  will  be 
handled  in  the  DSM  instead  of  complicating  the  decomposition. 


5.2  Design  Structure  Matrix  (DSM) 


5.2.1  General  DSM  Mechanics 

The  DSM  method  was  chosen  to  synthesize  the  current  architecture  because  of  the 
existing  modular  design  of  the  first  generation  NILM  concept  and  element  interconnections.  For 
a  quick  review  of  how  DSM  works,  consider  a  system  that  is  composed  of  two  elements  /  sub¬ 
systems  (or  activities/phases):  element  "A"  and  element  "B".  A  graph  may  be  developed  to  represent 
this  system  pictorially.  The  graph  is  constructed  by  allowing  a  vertex/node  on  the  graph  to  represent 
a  system  element  and  an  edge  joining  two  nodes  to  represent  the  relationship  between  two  system 
elements.  The  directionality  of  influence  from  one  element  to  another  is  captured  by  an  arrow  instead 
of  a  simple  link.  The  resultant  graph  is  called  a  directed  graph.  There  are  three  basic  building  blocks 
for  describing  the  relationship  amongst  system  elements;  parallel  (or  concurrent),  sequential  (or 
dependent)  and  coupled  (or  interdependent)  as  shown  in  Fig.  5-5. 


I  III  re  (  oiingiii  atiniis  tiiuf  (  lisirurtrii/r  :i  System 

^  Relationship 

Parallel 

.Sequential 

(  oiiplrd 

Grapii 

Repiesentalioii 

■CqI* 

[§] 

Figure  5-5:  DSM  relationship  building  blocks 


The  typical  matrix  representation  of  a  digraph  is  a  binary  (i.e.  a  matrix  populated  with  only 
zeros  and  ones)  square  (i.e.  a  matrix  with  equal  number  of  rows  and  columns)  matrix  with  m  rows 
and  columns,  and  n  non- zero  elements,  where  m  is  the  number  of  nodes  and  n  is  the  number  of  edges 
in  the  digraph.  The  matrix  layout  is  as  follows:  the  system  elements  names  are  placed  down  the  side 
of  the  matrix  as  row  headings  and  across  the  top  as  column  headings  in  the  same  order.  If  there  exists 
an  edge  from  node  i  to  node  j,  then  the  value  of  element  ij  (column  i,  row  j)  is  unity  (or  flagged  with 
a  mark  such  as  “X”  or  “•”).  Otherwise,  the  value  of  the  element  is  zero  (or  left  empty).  In  the  binary 
matrix  representation  of  a  system,  the  diagonal  elements  of  the  matrix  do  not  have  any  interpretation 
in  describing  the  system,  so  they  are  usually  either  left  empty  or  blacked  out  as  shown  in  Fig.  5-6 
[32]. 
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Three  Configiiratious 

that  Characterize  a  System 

Relationship 

Parallel 

Sequential 

Coupled 

DSM 

Representation 

A 

B 
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A 

■ 

A 

■ 

A 

■  x 

B 

■ 

B 

X  ■ 

X  ■ 

Figure  5-6:  Binary  representation 


5.2.2  DSM  with  expanded  connections 

Working  with  Oliver  de  Week,  Professor  of  Aeronauties  and  astronauties  and 
Engineering  Systems,  a  DSM  using  an  expanded  eonneetion/interfaee  definition  was  applied  to 
the  first  generation  NILM  system.  This  DSM  was  built  with  four  sub-cells  that  capture  the  four 
primary  types  of  connections  /  interfaces  common  in  electromechanical  systems.  The 
connections  between  the  elements  are  shown  as  the  off-diagonal  elements.  Those  connections 
include:  physical  connection,  mass  flow,  energy  flow  (includes  electrical,  mechanical,  or 
thermal  flow),  and  information.  Figure  5-7  presents  a  simple  example  of  how  to  employ  the 
expanded  DSM. 


A  sinnple  system 


Block  Diagram 


Design  Structure  Matrix 


Figure  5-7:  Expanded  Connections  DSM  Example  132] 


Component  A  physically  connects  to  B  which  in  turn  is  connected  to  C.  A  mass  flow 
occurs  from  B  to  C,  while  energy  is  supplied  from  A  to  B  and  C,  respectively.  Additionally,  A 
and  B  exchange  information  with  each  other.  Such  a  DSM  forms  the  basic  information  upon 
which  the  subsequent  analysis  builds. 


Each  of  the  level  2  parameters  was  entered  into  the  expanded  DSM  according  to  the 
inter-relationships  described  below.  The  goal  of  this  DSM  will  be  to  analyze  and  critique  the 
first  generation  architecture  and  suggest  an  improved  architecture.  Table  5-2  represents  an 
expanded  DSM  with  first  generation  NILM  architecture  entered  prior  to  partitioning. 
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Table  5-2:  DSM  with  expanded  Interconnections  /  Interfaces 


NILM  Baseline 
DSM 

Configuration 


Each  position  that  would  ordinarily  represent  a  standard  conneetion  or  interface  has  been 
sub-divided  into  four  positions.  Starting  in  the  upper  left  corner  and  working  elockwise  the 
positions  represent  Physieal  conneetion,  Mass  flow,  Information  flow,  and  Power  flow.  The 
NILM  system  in  its  current  form  does  not  have  any  Mass  flow  eonnections.  In  order  to  keep  the 
data  entry  and  subsequent  partitioning  as  simple  as  possible,  some  standards  for  data  entry  were 
observed. 

They  are  as  follows: 

•  If  energy  or  information  flow  depended  on  using  a  standard  interfaee  both  physical 
connection  and  flow  connection  was  entered  in  the  eorresponding  intersection. 


78 


•  Any  design  dependency  was  considered  a  two  way  dependency.  This  is  always  the  case 
for  physical  connections  hut  is  also  true  for  some  power  and  information  sub-systems. 
For  example,  the  PC  power  supply  depends  on  the  UPS  enclosure  for  a  connection  and 
the  size  of  the  UPS  enclosure  depends  on  PC  power  requirements. 

Effort  was  made  to  ensure  all  like  subsystems  are  made  dependent.  For  example,  all 
similar  power  supplies  are  made  dependent  to  one  source.  The  design  team  also  decided  that  the 
NILM  system  should  have  the  ability  to  draw  power  from  the  target  system  instead  of  needing  a 
separate  power  source.  For  example,  most  target  shipboard  systems  use  three  phase  440V 
power.  A  step-down  transformer  of  some  kind  could  be  used  to  deliver  1  lOV  single  phase  power 
and  eliminate  a  requirement  to  have  extension  cords  run  across  the  engineering  spaces.  In  order 
to  represent  this,  a  physical  dependency  was  established  prior  to  partitioning  between  the  target 
system  and  the  single  phase  power  used  to  power  NIFM  loads.  With  these  adjustments  made, 
the  matrix  was  partitioned. 

5.2.3  Partitioned  DSM Results 


Table  5-3:  Partitioned  DSM 
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Table  5-3  represents  the  partitioned  DSM.  The  partitioned  matrix  presents  a  reduced 
number  of  sub-systems  or  reduced  level  1  decomposition,  now  there  appear  to  be  three  sub¬ 
systems.  The  main  subsystem  in  the  upper-right  has  combined  the  PC,  Display,  and  UPS  sub¬ 
systems,  followed  by  the  target  or  monitored  system  and  finally  the  new  NEMA  enclosure 


5.3  Implementing  Results 

5.3.1  COTS  Hardware 

At  this  point  the  team  began  investigating  how  they  could  move  forward  with  selecting 
hardware  to  fulfill  the  needs  of  each  sub-system.  COTS  systems  have  been  successfully 
employed  to  handle  the  processing  and  display  functions  aboard  ship.  The  team  met  with  the 
end  users  and  engineering  leadership  to  determine  specific  needs.  At  the  same  time,  diagnostic 
software  was  being  developed  that  would  enable  near  real  time  monitoring  of  target  systems. 
The  resulting  needs  of  the  new  COTS  hardware  component  were: 

•  Able  to  run  Linux  and  Windows  environment 

•  Small  footprint  in  the  engineering  spaces 

•  Out  of  the  way 

•  Water  resistant 

•  Ethernet  compatibility 

•  Touch  screen 

•  Easy  to  mount  /  and  remove 

The  design  team  felt  that  a  tablet  PC  could  accomplish  the  task  of  the  first  partition.  A 
touch-screen  display  and  stylus  could  be  used  as  the  primary  user  interface.  The  team  decided 
on  a  Lenovo  X61  tablet  PC  because  of  its  proven  design  and  cost.  Additionally,  accessories 
were  available  that  would  assist  in  protecting  the  device  in  arduous  environments  like  those 
found  aboard  ship. 


Figure  5-8:  Lenovo  X61  Tablet  PC  and  Tablet  Carrying  Sleeve  (bttp://www.leno vo.com) 
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To  protect  the  tablet,  the  team  decided  on  modifying  a  case  similar  to  the  one  pictured  to 
the  right  in  Fig.  5-8.  The  PC  is  secured  inside  the  sleeve  with  Velcro  flaps  that  can  protect  the 
unit  from  minor  debris  and  indirect  water  spray.  The  clamping  mechanism  used  to  secure  the 
straps  were  very  durable  so  the  team  decided  to  simply  remove  these  from  the  straps  and  bolt 
them  to  standard  mounting  brackets  typically  used  to  secure  light  equipment  to  engineering 
bulkheads.  This  tablet  /  PC  sleeve  combination  met  the  fulfilled  the  requirements  of  the  COTS 
partition. 


5.3.2  New  NEMA  Enclosure 

The  most  significant  changes  suggested  by  the  DSM  involved  the  NEMA  enclosure.  The 
advances  in  COTS  hardware  and  need  for  increased  capability  drove  some  major  changes  in  the 
NEMA  enclosure.  The  vision  for  the  second  generation  NILM  included  the  ability  to  use 
streaming  data  for  near  real  time  diagnostics.  It  needed  to  be  easier  to  install  and  initialize  than 
past  systems  for  commercial  appeal.  The  resulting  needs  of  the  NEMA  Enclosure  were: 

•  Minimal  footprint 

•  Ability  to  stream  up  to  six  channels  of  data 

•  Simplification  of  transducer  scaling  resistor  sizing 

•  Ethernet  ready 

•  Ability  to  power  itself  off  typical  marine  440V  system 

The  result  was  the  new  second  generation  NILM  NEMA  enclosure  pictured  in  Eig.  5-9. 


NEMA 
Power  - 
supply 


Universal 

Transformer 


I _ ^ _ 

Figure  5-9:  Second  Generation  NILM  NEMA  Enclosnre 
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Lab  Jack 
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The  universal  transformer  can  be  configured  to  step  208  -  480  Volts  down  to  a  100  -  125 
Voltage  range  that  can  be  used  by  the  NEMA  power  supply.  The  new  enclosure  can  support 
monitoring  of  both  delta  and  Wye  configurations. 

In  the  past,  current  transducer  resistor  sizing  was  very  tedious.  A  power  meter  was  used 
to  determine  a  range  of  current  the  system  to  be  monitored  drew  from  the  load  center.  Then  a 
calculation  was  made  to  determine  the  sizing  resistors  that  were  to  be  installed  across  the 
terminal  block  in  Fig.  5-1.  With  these  sizing  resistors  organized  in  a  DIP  (Dual-Inline  Package), 
an  installer  can  quickly  select  a  safe  value,  run  an  initial  test  and  immediately  set  the  current 
transducer  sizing  resistors.  The  new  software  developed  by  LCDR  Ethan  Proper  and  PhD 
candidate  Jim  Parish  can  scale  both  current  and  voltage  transducers. 

One  significant  detail  that  fell  out  of  the  DSM  was  that  the  Analog  to  Digital  (A  /  D)  converter 
should  be  moved  outside  the  PC.  The  team  began  looking  for  a  device  that  could  be  installed 
into  the  NEMA  enclosure.  Additionally,  a  high  data-rate  or  one  that  could  perform  the  A/D 


Figure  5-10:  Lab  Jack  Model  UE9 

function  at  a  rate  of  at  least  48kbits  /  second  was  desired  (six  channels  sampled  at  a  rate  of  8kbits 
/  second)  to  make  streaming  data  possible.  The  Lab  Jack  UE9,  USB/Ethernet  based 
multifunction  data  acquisition  and  control  device  shown  in  Fig. 5- 10  was  chosen  to  handle  the  A  / 
D  function.  While  this  product  seems  more  complex  than  needed,  it  enables  data  streaming. 
This  made  near  real  time  interactive  diagnostics  possible. 


5.3.3  Target  System 

The  target  system  or  system  to  be  monitored  will  now  also  supply  power  to  most  of  the 
NILM  equipment.  Voltage  connections  and  current  transducer  installation  will  generally  be 
done  the  same  as  with  the  first  generation  NILM.  One  change  involves  configuring  the  voltage 
block  seen  in  Fig.  5-9  to  support  host  system  power.  Two  basic  issues  must  be  considered  with 
respect  to  the  voltage  block.  The  first  is  whether  a  neutral  is  used  (single  phase  or  Wye 
configurations)  and  the  second  is  choosing  the  phase  of  host  power  to  be  used. 
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5.4  Completed  Installation 


The  new  NILM  hardware  and  software  was  developed  and  tested  in  LEES  laboratory 
using  a  simple  motor  and  dynamometer  setup.  Onee  the  team  was  eomfortable  with  the  second 
generation  NILM,  it  was  installed  aboard  ship.  The  installation  in  Eig.  5-3  was  removed  and  the 
second  generation  NILM  installed.  Eigure  5-11  was  taken  aboard  Escanaba  shortly  after 
installation  and  testing.  Note  how  much  cleaner  the  installation  looks. 


5-11:  Second  Generation  NILM  Install  aboard  CGC  ESCANABA 
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5.4.1  General  Installation 


Installation  of  a  second  generation  NILM  set  is  quite  different  from  its  predecessor  and 
less  complicated  for  the  most  part.  However,  because  of  the  additional  capabilities  of  the  new 
system,  target  system  operational  understanding  becomes  more  essential.  The  NILM’s  impact 
on  the  system  must  be  taken  into  account.  This  section  is  not  meant  to  be  a  complete  installation 
guide  but  is  meant  to  provide  a  general  outline  for  installation. 


5.4. 1.1  Survey  Target  System  Electrical  interface. 

The  first  step  to  installing  a  NILM  system  on  a  particular  piece  of  equipment  is  to  obtain 
the  load  panel  electrical  line  drawings.  This  is  important  for  a  number  of  reasons.  The  first  of 
which  is  to  locate  a  position  downstream  of  the  last  possible  protection  device.  Figure  5-12  is 
the  inside  of  the  CHT  load  panel  after  NILM  install.  The  area  circled  in  red  represents  the  best 
place  sample  voltage  and  install  current  transducers. 


Figure  5-12:  Open  CHT  Load  Panel 


A  power  meter  should  be  used  during  the  survey  process  to  observe  a  complete 
operational  cycle  and  maximum  observed  current  noted.  LEM  current  transducers  (CTs)  use 
model  numbers  like  LA-150P,  the  numerical  value  in  the  model  number  is  the  maximum  current 
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the  CT  is  designed  to  operate  with.  Once  the  proper  CT(s)  is  chosen,  the  target  system  must  be 
studied  to  determine  the  optimal  phase(s)  to  monitor.  The  equipment  line  diagram  will  come  in 
handy  here  as  well.  For  example  Fig.  5-13  presents  a  power  reading  taken  on  phase  A  with  the 
system  in  standby.  Note  that  phase  A  shows  a  200  Watt  load.  In  some  systems,  200  Watts  is 
enough  to  cause  a  software  diagnostic  package  hardship. 


5-13:  Auxiliary  Load  on  CHT  System  Phase  A 


In  the  case  of  the  CHT  system,  the  best  channel  to  monitor  would  be  phase  C.  The  next 
thing  to  consider  is  how  best  to  access  the  load  panel.  Ships  usually  do  not  want  extra  holes 
drilled  in  equipment  enclosures  unless  absolutely  necessary.  Any  installation  must  maintain  the 
integrity  of  the  enclosure.  Meaning,  accesses  must  be  semi-water  tight. 


Figure  5-14:  Accessing  the  Load  Center 


85 


Once  the  CT  and  voltage  connections  have  been  made  and  properly  passed  through  the 
target  system  enclosure,  the  NILM  voltage  block  must  be  configured  depending  on  the  type  of 
voltage  sampled.  NILM  boxes  are  manufaetured  three  phase  Wye  ready.  However,  most  ships 
are  three  phase  delta.  The  best  plaee  to  make  a  delta  to  Wye  conversion  is  at  the  voltage  block 
circled  in  red. 


Figure  5-15:  Second  Generation  NEMA  Voltage  Terminal  Block 


Mount  the  Computer  in  a  loeation  out  of  the  way  at  eye  level.  It  is  usually  best  to  survey 
the  system  end  users  to  determine  best  plaeement.  The  NILM  NEMA  enclosure  should  be 
mounted  as  close  to  the  load  eenter  as  possible  in  a  manner  not  to  interfere  with  surrounding 
equipment  and  at  least  12  inches  off  the  deek. 


5.5  Conclusions 

As  NILM  moves  forward,  it  is  important  for  it  to  evolve  into  a  system  that  ean  support  a 
wider  eustomer  base  and  yet  maintain  the  flexibility  for  employment  in  both  shipboard  and  shore 
based  applieations.  The  new  NILM  arehiteeture  offers  an  ineredible  leap  in  capability  and  value 
to  customers.  The  arehiteeture  should  be  revisited  every  two  years  to  ineorporate  the  benefits  of 
new  COTS  technology  as  well  as  host  ship  systems. 


86 


6  Conclusions  and  Recommendations 


6.1  Single  Machine  Monitoring 

NILM  has  shown  the  ability  to  monitor  a  single  induction  machine  using  current 
harmonic  quantities.  However,  there  are  both  equipment  orientation  and  physical  design 
characteristics  that  can  enhance  this  process.  It  was  found  that  a  simply  loaded  machine  is  more 
easily  monitored  in  a  horizontal  direction  vise  a  vertical  one.  Machine  loading  can  potentially 
dampen  the  spectral  quantities  as  well.  These  are  considerations  when  monitoring  both  single 
machines  and  aggregate  load  centers. 

A  single  machine  diagnostic  package  should  be  developed.  A  suggested  architecture  for 
such  a  system  is  presented  in  Chapter  3.  The  basic  physical  building  blocks  for  such  a 
diagnostic  system  would  be  very  similar  to  the  ones  presented  in  Chapter  5.  A  simple  system  or 
induction  machine  should  be  used  to  initially  test  the  system. 


6.2  Multiple  Machine  Monitoring 

NILM  has  the  ability  to  monitor  multiple  induction  machines  using  an  aggregate  sensor. 
The  test  outlined  in  chapter  4  of  this  thesis  should  be  replicated  using  two  induction  machines 
with  differing  key  parameters  and  increased  complexity.  A  system  like  the  automated  variable 
speed  drive  main  engineroom  ventilation  system  used  aboard  the  Marine  Protector  Class  Cutter 
should  be  investigated.  This  would  present  the  NILM  with  a  two  induction  machines  and  solid 
state  variable  speed  drive. 


6.3  Second  Generation  NILM  Architecture 

The  architecture  presented  in  this  thesis  was  made  possible  primarily  because  of  COTS 
hardware  affordability.  Five  years  ago  the  team  would  not  have  felt  it  was  economically  feasible 
to  put  a  touch  screen  tablet  PC  on  a  Military  vessel.  Additionally,  new  product  architectures  like 
the  Lab  Jack  (UE9)  multifunction  data  acquisition  and  control  module  that  enabled  near  real  time 
data  streaming  as  well  as  analog  to  digital  functions.  The  current  architecture  should  be 
reviewed  in  the  next  two  years  to  take  advantage  of  affordable  technologies  as  well  as  ship 
technology  advancements.  For  example,  the  ability  to  remotely  transmit  diagnostic  information 
should  be  investigated. 
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